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Abstract 
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Tbs US6S of radar observations of planets and very—long'—baseline 
intorferoniGtric observations of extragalactic objects to test theories of 
gravitation are described in detail with special emphasis on sources 
of error. The accuracy achievable in these tests with data already 
obtained, can be summarized in terms of 
standard notation as follows : 

Retardation of signal propagation 
(radar) : 

Deflection of radio waves (inter- 
ferometry) : 

Advance of planetary perihelia 
(radar) : 

Gravitational quadrupole moment of 
sun (radar) ; 

Time variation of gravitational 
constant (radar) : 


a (ill, 

0.04 

a (ill) 

0.03 

,(2+2Y-e) 

0.04 

a(J2> 

• 

1.5 X 10"^ 

a(§) 

-1 

= 10 yr 


The analyses completed to date have yielded no significant disagree- 
ment with the predictions of general relativity. If radar 
and radio observations are pursued in the manners herein proposed, the 
uncertainties in these experiments could be reliably reduced to 


.(ill) 

< 

** 

0.0Q3 (radar) 


< 

0.003 (radio interferometry) 



0.01 (radar) 


< 

3 X 10 ^ (radar) 

“(|) 

< 

2 X 10 ^^yr ^ (radar) 


biy the early 1980's. In addition, a significant radar test could be 
made in this time period of the Principle of Equivalence, in regard 
to the relative contributions of gravitational binding energy to 
inertial and gravitational masses. Combination of radar and space- 
craft ranging data would yield improvements in almost all of these tests 
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Foreward 


■ ■ We are concerned in this report with the evaluation 

of the state of the art in testing theories of gravitation 
using two radio-astronomy methods: radar time-delay measure- 

ments and very-long-baseline interferometry (VLBI) ray-deflection 
measurements. For convenience, we have therefore divided the 
report into two parts. In Part I we discuss the radar 
method, but give relatively little emphasis to a description 
of the basic measurement technique since it is already well 
documented in the literature. In Part II, by contrast, we 
describe in some detail the VLBI technique as applied to testing 
theories of gravitation, since this method is newer and con- 
sequently less well known than the radar method. In both 
parts of the report, we conclude with a consideration of 
possible improvements in the accuracy of the relevant tests of 
gravitation. 
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Part I 

Radar Tests of Relativity 

I, Introduction 

Radar tests of the theory of general relativity, or, more 
generally, of theories of gravitation have been in progress 
for more than a decade and several useful results have already 
been obtained. Reviews of this work, containing references 
to the original publications, can be found in several 
articles (Shapiro, 1972 and 1973) . 

Interplanetary radar measurements which form the basis 

of these tests primarily concern the round-trip time delays 

and the Doppler shifts of the echoes. The time delays have 

fractional uncertainties far less than those for the Doppler 

shifts and are the principal measurements underlying the 

tests of theories of gravitation. Relativistic effects, on 

general grounds, can be expected to appear when the fractional 

( 2 / 2 . 

uncertainties in the delay measurements dip below vv /c j 

lO”^, with V a typical orbital velocity for a planet, 

—7 

or below about 10 for ray paths that nearly graze the solar limb. 

Radar system sensitivity has been increasing, on average, by 
about 5.5 db/yr since 1946 and, at present, fractional delay- 
measurement uncertainties for some favorable planetary con- 
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figurations are as low as 2 parts in 10^®. On the horizon is 
the upgraded Arecibo radar facility which, after it begins 
routine operation in mid 1975, should extend such ac- 
curacies to most configurations of the inner planets as well 
as increasing the accuracies for certain configurations. The limitations 
in these latter circumstances will be set by considerations 
other than the basic radar system sensitivity. 

The above statements indicate the general feasibility of 
performing precision tests of theories of gravitation with 
radar measurements. However, one must be cautious: In ad- 

dition to information about relativistic effects on the signal 
propagation and on the motion of the observed body and the 
observing platform , there exist other nonrelativistic in- 
fluences on the delays. These can be either of a random, un- 
predictable nature or of a type that has* a signature highly 
correlated with the sought- for relativistic effects. In 
both cases, the net result is to degrade the accuracy of the 
tests of gravitation theories and, therefore, such other in- 
fluences will be considered "corrupting" for the purposes of 
this report, 


We divide the remainder of this part into three main 
sections. In the first, we discuss the various radar tests 
that appear feasible and include approximate magnitudes of the 
relativistic effects involved. In the second, we consider in 
a rather encyclopedic fashion, effects that might corrupt the 
tests and the means for the elimination or reduction of their 
adverse influences. The final section is devoted to a brief 
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summary of the present status of the relativity tests and an 
indication of the possible improvements that could be anti- 
cipated in the next decade. 

II. Relativistic Effects on Radar Measurements 

We discuss relativistic effects primarily, but not ex- 
clusively, in the context of the Schwarzschild solution. With 
one exception, we consider the relativistic effects of the 
sun to stem from a non-rotating, spherically-symmetric mass 
distribution. The relativistic effects stemming from the non- 
zero masses of the planets can be ignored. This last con- 
clusion follows from the formal solution to the post-Newtonian 
(weak field) equations for the solar system, derived many times 
by many authors over the past half century (see, for examples, 
Tausner, 1966 and Weinberg, X912) , 

Our direct concern is with the relativistic effects caused 
by the sun on the propagation or radio waves and on the motions 
of the inner planets , Mercury through Mars . The orbits and 
masses of the outer planets may be assumed known in that the 
uncertainties in these quantities, determined separately from 
the direct observations of the outer planets , cause uncertain- 
ties in our interpretation of the inner-planet observations 
that with few exceptions (see Section III) are too small to af- 
fect the relativity tests presently of interest. 

1, Parameterization of the Schwarzschild Metric 

Eddington (1960), and later Robertson (1962), sought 
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to generalize the space-time metric that represents the 
(exterior) Schwarzschild solution for the sun partly to 
formalize comparisons with experimental results. They chose 
a particular parameterization which is analogous to the 
parameterized multipole expansion often used for an unknown 
charge distribution. For the Schwarzschild solution ex- 
pressed in so-called isotropic coordinates, the Eddington- 
Robertson form of the metric is 

IT r 

ds^ = [l-2a(^) + 26 + ...Ic^dt^ 

’-[1+2y(^) + . . . ] [dx^+dy^+dz^] , (2.1) 

where ds represents the element of arclength in the 
four-dimensional space-time; dx, dy, and dz are the dif- 
ferentials of the space coordinates; dt is the differential 
of time; c is the speed of light; and a, 6f and y the para- 
meters.* 

A more complete, and very elegant, parameterization 
of metric theories of gravitation in the post-Newtonian 
approximation has been developed over the past few years by 
Nordtvedt, Thorne, and Will (see Will, 1973 for a thorough 
discussion). However, the parameterization given in Equation (2.1) 


In view of the excellent agreement of planetary orbits with 
Newtonian theory, for example, one can conclude without further 
ado that a is equal to unity. 
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will, for the most part, be adequate for our purposes. 

Since the parameters y and B in this parameterisation 
affect both signal propagation and planetary motions, we 
will discuss each effect in turn,- 

2 . Retardation of Signa l Propagation 

The Eddington-Robertson metric, as does the Schwarz- 
child solution itself, predicts a retardation of, or delay 
in, signal propagation when rays pass near a massive body, 
in particular, as was first shown a decade ago (Shapiro, 1964), 
the "excess" delay, Ai^, is predicted to obey: 


2(l+y)ro 

Atr= c tg+rp-R ' 


( 2 . 2 ) 


where r^ 5 (GM^/c^) = 1.5 1cm is the gravitational radius 

of the Inn, G the gravitational constant: the sun's mass: 

o the speed of light: r^ the heliocentric distance of the 
earth: r the heliocentric distance of the target planet: and 
R the earth-planet distance. The maximum predicted effect (with 
y=l) is approximately 250 pseo for a signal path that grazes 
the solar limb. But, as intimated above , the "observable" 
effect is, in effect, reduced by other nonrelativistic in- 
fluences on the round-trip time delays. 
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Note that the expression {2.2), accurate to first order 
in r^, is independent of 3. A dependence of the time delay 
on 3 enters to this order only through the relativistic ef- 
fects on planetary orbits. The direct effects of solar gravity 
on the time delays are sensitive only to y (Shapiro, 1966), 

3. Deflection of Signals 

Interferometric radar observations of a planet can 


be employed to measure the predicted gravitational de- 
flection, Ti, of the radar signals by the sun (Shapiro, 1967)^ 


(1+Y)2^ 


n = 


° tan(|) , 


(2.3 ) 


where 0 is the planet-sun-earth angles Because of the finite 
planet-earth separation, the maximum deflection of 1.75 arcsec 
(Y=1) is not attained. Thus, for Venus the target planet and 
with the signals grazing the limb of the sun, n - 0.73 arcsec 


(Y= 1 ) . This experiment has never been performed because higher 
accuracy appears to be obtainable from the corresponding radio 
interferometric observations of cosmic radio sources, as dis- 
cussed in Part II of this report. 

4. Advance of Planetary Perihelia 

As has been known since the inception of the theory 
oyer a half-century ago, general relativity and its Ed- 
dington-Robertson generalization predict a secular advance of 
planetary perihelia, given by 

6 TT 1C 

Am - ( ^•*'^T~' J-) — S. rad/rev, (2.4) 

3 p 

where £ is the semilatus rectum of the orbit. Clearly the 
advance is sensitive to both y and 3 with the sensitivity to 
y being twofold greater. The effect of this excess ad- 
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vance on the measurement of time delays is boianded, very 
crudely, for radar observations of Mercury by 

Aiadv 2eaAo) - 160 ysec/yr (2.5) 

where e is the eccentricity and a the semimajor axis ox 
Mercury's orbit. The value- for Aw is given by Equation (2.4) 
but scaled to the total advance per year, about 0.4 arcsec; 
both Y and 6 were set to unity in accord with the predictions 
of general relativity. The factor of tv/o in Equation (2.5) 
accounts for the round-trip aspect of the delays. 

There are other, short-period relativistic terms in the 
description of planetary orbital motions. These, for 
Mercury, have amplitudes of up to about 40 Usee. But their 
"observability" is doubtless much diminished by correlations 
with the orbital elements, or initial conditions, of planetary 
motion which perforce must also be estimated from the same 
time-delay data since no other data are of sufficient accuracy. 

A similar comment applies to the relativistic deviations pre- 
dicted from Kepler's third law for a suite of planetary orbits. 
As far as we are aware, no careful study has yet been carried 
out to isolate, in an operational manner, the contributions each 
of these orbital effects could make to the determinations of 
Y and 6, and, in particular, to their separation from the 
^y^^^ical effects of the solar gravitational quadrupole moment 
(see Section Ill.iii.l) . in all covariance studies so far 
carried out all such effects are lumped together. 

5. Possible Time Variation of the Gravitational Constant 
According to both Newtonian precepts and Einstein's 
theory of general relativity (and its Eddington-Robertson 
geheralization) , the constant of gravitation is indeed a 
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universal constant. On the basis of other theories of gravi- 
tation, such as the Brans-Dicke scalar-tensor theory 
(see, for discussions, Weinberg, 1972 and Misner, Thorne, 

t . 

and Wheeler, 1973) , the constant of gravitation should be 
weakening slowly with time as the universe expands. This 
weakening would be observable through measurements that 
depended on atomic constants. In other words, on the basis 
of such theories, a "gravitational" clock would appear to slow 
down as measured by an atomic clock. In particular, if the 
mean motion, n, of a planet were measured with an atomic 
clock, it' would be predicted to decrease in a secular fashion 
according to the formula 



( 2 . 6 ) 


where a dot indicates differentiation with respect to time. 
This effect would be the more noticeable the larger the 
mean motion. Thus, of the planets. Mercury is the most 
useful gravitational clock for the detection of any such 
time variation of G. The moon, although it has approximately 
a threefold higher mean motion than Mercury, suffers from 
being too close to the earth with the consequence that tidal 
interactions are strong and difficult to separate from the 
effects of any possible change in G. 

A change in G causes a change in the earth-planet time 
delay, due to a change, AL, in the mean longitude of the 
planet, bounded by (Shapiro et al., 1971a): 
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At, 


- 2aAL = 2an^T^ , 


(2.7) 


since 


T. , 

AL = / n t dt = ^ nT 
0 ^ 


2 


G „2 
^ G ^ 


( 2 . 8 ) 


Here T is the interval over which the delay observations 

-10 

have extended. For Mercury, with (G/G) -lx 10 per 
year, we obtain 

Ata - ysec (2.9 ) 

with T in years. Thus At^ - 25 ysec after 5 years. 

Note that there have been attempts to estimate G/G from 
studies of lunar motion, notwithstanding the difficult problems 
posed by tidal interactions (see, for example, Slade 1971). 

In the most recent attempt. Van Flandern (1975 ) claims to have 
detected a significant change: (G/G) = (-7.5 +2.7) x 10 per 

year, based partly on an analysis of 18 years of lunar occul- 
t^tion data which allow him to estimate the acceleration m mean 
longitude shown in Equation (2.8) . However, 
for the moon, in 18 years, the change in AL would be only 
about 0.4 arcsec for this value of (G/G). 
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Again we must: emphasize that the whole effect here is 
by no means ’’observable". First, we note that a parabola, 
such as given by Equation (2.9), can be fit by a straight 
line with the deviation nowhere exceeding one eighth of the 
total "growth" of the parabola. Thus, unless measurement 
errors were under, say, 3 usee, we would have difficulty, 
based on' only 5 years of earth-Mercury time-delay observa- 
tions, in distinguishing a change in G of 1 part in 10^^ 
per year from a modification in our estimate of Mercury's 
mean motion. [Recall that the slope of the best-fit straight 
line to the parabola of Equation (2.9) would be interpreted, 
in this oversimplified analysis, as a correction to the 
estimate of Mercury's mean motion.] Second, we mention 
that the actual estimate of G would involve a simultaneous 
solution for all relevant parameters with the consequence that 
the "observable" effect would actually be reduced by an amount 
even larger than the factor of eight reduction entailed by the 
effective "masking" of a parabola that can be accomplished 
by a straight line. 
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The abqve arguments apply as well 
to the interpretation of changes in the moon's mean longitude 
that might be due to a variation in the gravitational con- 
stant (see Section IV. 4) . In view of these and other 
aspects of the analysis, we consider that Van Flandern has 
underestimated the uncertainty of his result which, in our 
opinion, is not significantly different from zero. 

Finally, we remark that in certain theories of gravi- 
tation, spatial variations of the gravitational constant 
are predicted. It does not appear feasible to use inter- 
planetary radar measurements to detect or to set useful 
bounds on such possible variations. However, a detailed 
analysis of this problem has not been carried out. 

6. Possible Violation of the Principle of Equivalence 
A cornerstone of Einstein's theory of general 
relativity is the Principle of Equivalence. In its so-called 
wfeak form, this principle states that the ratio of the 
gravitational to the inertial mass of a body is independent 
of its composition and size. This principle has been tested 
in the laboratory with ever-increasing accuracy over the past 
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three centuries, starting with Newton's demonstration which 
established the validity of the principle to about 1 part in 

3 

10 . The most accurate test was carried out by Braginski and Panov 

(1971) who concluded that the principle held to at least a 
12 

few parts in 10 . Although these laboratory tests demonstrate, 

for example, that the binding energy of the electrons in 
atoms and of the nucleons in the nucleus contribute equally 
to inertial and gravitational masses, these tests fail utter- 
ly to shed light on whether the gravitational binding energy 
contributes equally to inertial and gravitational masses; 

For a laboratory-sized object, the gravitational binding 

23 

energy constitutes no more than about 1 part in 10 of the 
total mass, far below the sensitivity of the laboratory tests. 
Since this fractional contribution varies approximately as the 
square of the body's dimension, one needs planetary-sized 
bodies to perform a meaningful test for the gravitational 
binding energy. Let us define the ratio of the gravitational 
to inertial mass such that 



where A, by assvimption, accounts for any possible deviation 
from the Principle of Equivalence due to an unequal contribu- 
tion of the gravitational binding energy to and M^. For 
the sun, then, we'd expect A to be bounded by O (10 cuid. 
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—8 —9 

for Jupiter and the earth, by about 0(10 ) and 0(10 ), res- 

pectively. 

A two-body orbiting system is not useful for a test 
of the principle of equivalence since the implied violation 
of Kepler’s third law, which relates the period and semi- 
major axis of such a system, is unobservable due to the 
lack of an independent determination of the relevant masses . 
However, three or more mutually orbiting bodies do provide 
a possible test. 

Renewed interest in such a test was kindled by Nordt- 
vedt (1968) who showed that this weak Principle of Equi- 
valence was violated for massive bodies in the Brans-Dicke 
theory, a consequence not previously recognized. More speci- 
fically , Nordtvedt showed that a violation would introduce 
a variation in the radial separation of two of the bodies 
(say, the earth and the moon) of a three-body system (say, 
the earth-moon- sun system) proportional to: 

48 - Y - 3 - + 1^2 - - j?2- (2.11) 

We, omit the exact definitions of the parameters a 2 » 

and ^ 2 ' they are not germane to our discussion; see 
Will (1973) for a detailed description. We remark here 
only that, for general relativity, the terms 
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on the right side of Equation (2.11) sum to zero whereas 
this sum does not vanish under the Brans-Dicke theory. 

Considerable attention has been given to the possible 
detection of a nonzero value for A from analysis of lunar 
laser observations, the relevant three-body system being 
the earth, the moon, and the sun. In the radar context, 
we have examined the possibility for detection of a 
violation with respect to the four-body system: earth-Mars- 

Jupiter-sun. In particular, we have considered using earth-Mars 
time-delay observations; a very detailed, but still incomplete, 
study of this possibility wasmade at M.I.T. by Sherman (1973). 

similar observations of Venus , especially , and of Mercury will also be 
very efficacious in strengthening this test. 

For a violation of the weak Principle of Equivalence of 
the order of the ratio of the gravitational self energy to the total 
energy of a body, it can be shewn that the effect on the earth-Mars time- 
delay neasurenents, for example, would be of the order of 10 usee (Sherman, 
1973) . But masking will reduce its observability. We note that the effect 
of a nonzero A will be in large part periodic, with periods approximately 

those of the earth and target planet. 

7, Gravitational ’’Itedshift" 

The gravitational "redshift", the effect of the local gravita- 
tional potenUal on the apparent frequency of a light wave vMch causes a 
spectral line generated at the sun to appear red-shifted when observed on 
earth, does not seem aitenable to direct. measurenent in interplanetary 
radar experiments. However, the concomitant predicted ef- 
fect on clock rates is susceptible of detection . Because 
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of the eccentricity of the earth's orbit, an atomic clock 
on earth, relative to that of a distant observer, is pre- 
dicted to undergo yearly variations of amplitude about 
1.5 msec. This effect, relative to the situation in its 
absence, causes the time-delay measurements to be referred 
to a different epoch. If the reference epoch is in error, 
however, the calculated delay — in the absence of any other 
errors or inexact parameter values — would disagree with 
the observed delay by the product of the epoch, error and the 
time rate of change of the delay. For interplanetary ob- 
servations, the delay rate can be as high as 2 (v/c) - 200 ysec/sec 
Hence an epoch error of 1.5 msec would manifest itself in delay 
residuals of maximum magnitude of about 0.3 ysec. Such an ef- 
fect is too small to be reliably detected with useful accuracy at 
present, especially when masking is considered, and so will 
not be discussed further in this report. 

8. Lense-Thirring Effect 

Lense and Thirring (1918) predicted on the basis of 
the then new theory of general relativity that a rotating 
mass tended to rotate a "nearby" inertial frame in the same 
direction. Incorporation of the sun's rotation, therefore, 
produces a modification of the Schwarzschild solution which 
has consequences for both the propagation of radio signals 
and the orbits of planets. In addition, the spin axis of 

i 

these orbiting bodies would undergo, additional precession 
due to the Lense-Thirring effect. None of these consequences 
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of the theory of general relativity, or its generalizaf , 
so far been observed. 

A program undertaken at Stanford University, with NASA 
sponsorship, has for over a deoade been ooncerned with the 
development of a superconducting gyroscope to orbit the earth 
to detect the so-called "geodesic precession" and the much smaller lense- 
Thirring precession which, in the terrestrial context, aaounts to only seme 
hundredths of an arcsecond per year (see Misner, ttome, and Wheeler, 1973, 
and references cited therein) . Any actual measurements are still seme years 

in the future. 

■ A preliminary study was undertaken at H.I.T. (Miller, 1971) to ascertain 
the possible observability of the Lense-Ihirring effect on interplanetary time- 
delay measurements. The results, as expected, were not encouraging. The effect 
on the earth-^tercury time^elays, due to the effect on the orbit of Meroory, 
is of the order of 1 nsec per orbital revolution. Worse, the 
lack of precise knowledge of the sun's angular momentum vec- 
tor makes the effect even more difficult to detect. We 
therefore end discussion of this effect for the purposes of 

this report. 
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III. Effects Corrupting Relativity Tests 

Here we discuss the various effects on the radar time- 
delay measurements that, for the purpose of testing theories 
of gravitation, may be considered as corrupting influences, 
or as "noise". In contrast- to the discussion of relativistic 
effects, we treat these corrupting effects in approximately 
the reverse order of their magnitude, reserving the most 
serious — planetary topography — for last. 

1. Earth Rotation 

Since all observations are made from the earth's 

surface, the motion about the center of mass influences the 
tiine-delay observations. The earth's deviations from solid- 
body rotation with a constant angular velocity about a fixed 
axis in inertial space are conventionally described by the 
pirecession, nutation, polar motion, variations in UT.l, and 
sdlid-body tides of the earth. These deviations are currently 
kriown with sufficient accuracy that, for purposes of analyzing 
a i decade of interplanetary radar measurements, the uncertainty 
iri the position of a point on the earth's surface with res- 

I 

pect to an inertial frame is of the order of 5 meters or 
approximately 0.03 ysec in equivalent two-way light time. 



-25- 


Such an uncertainty is of no consequence for the interpretation • 
of existing radar data and will pose no severe problem for the 
interpretation of data to be obtained in the future. However, 
with the development of very-long-baseline interferometry 
(VLBI) , the uncertainty in the changes in the earth's position 
with respect to an inertial frame formed by distant extragalac- 
tic radio sources will be reduced over the next few years to 
the decimeter level (~^ 0.002 ysec in equivalent light time) and 
so will be an utterly negligible source of corruption for radar 
tests of relativity. 

2, Site Locations 

In addition to the determination of the earth's 
orientation in inertial space, one must also determine the 
location of the radar site on the earth. For some of the 
delay measurements already made, as well as for most of the 
future measurements , conventional geodetic surveys do not 
provide sufficient accuracy. Of course, one can add the site 
coordinates to the parameter list, as we have in fact done, 
and estimate these along with the other solar-system parameters 
of interest; the correlations with the estimates of the para- 
meters describing the relativistic effects are very low. 

If for no Other reason than as an independent check, how- 
ever, it is desirable to determine the site locations 
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by other means. The VLBI technique provides such a means; 
moreover, through its use, the radar sites can be located 

i , 

with respect to the same inertial frame in which the earth's 
orientation is described. We have already conducted a series 
of VLBI experiments between Golds tone and Haystack (Shapiro 

j 

et al . , 1974) which serves to establish their relative loca- 
tion to within 1 m*. We have also submitted a proposal to do 
a similar, but for technical reasons somewhat less accurate, 

VLBI experiment between Haystack and Arecibo. The results will allow the 
relative locations of all three sites to be determined to within a few meters 
which should provide a powerful check on the estimates of the site locations 
made from an analysis of the radar data from these three sites. 

We should also note that through VLBI the orientation 
of the planetary system itself can be determined with respect 
to the same inertial frame used to specify the earth's 


*From VLBI observations of extragalactic sources, the baseline 
of the interferometer is determined, with respect to the 
inertial frame formed by these sources, to within a parallel 
displacement. Lack of parallax leads to this freedom in 
the specification. 
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orientation (Counselman et al . , 1972)*. The method proposed to 
accomplish this determination involves differential measurements 
of extragalactic radio sources and transmitters aboard space- 
qraft that are either flying by, orbiting, or landed on anothc>r 
planet. By observing simultaneously, or near simultaneously, 
with a very-long-baseline interferometer, both the signals from 
the spacecraft and those from an extragalactic radio source, 
when the latter is nearly in the same direction from the earth 
as the planet, high accuracy can be obtained. Recently our VLBI 
group (see Part II) has used this technique to determine the 
angular separation between two extragalactic radio sources with 
an uncertainty of only about O'.'OOl; the same accuracy could be 
achievable in the determination of the orientation of the planets 
with respect to such extragalactic sources. The main problem is 
the- lack of suitable opportunities. 

3, Propagation Medium 

The medium through which the radar signals pass ob- 
viously affects the values of the measured time delays. 

For purposes of discussion, we separate the medium into 


*Because of the rather weak dynamical coupling between the 
earth's spin angular momentum vector and the angular momentum 
vector of the planetary system, the determination of the 
orientation of one of these vectors with respect to a given 
frame doesn't of itself determine the orientation of the 
other with respect to that frame. 
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three parts: the earth's atmosphere and ionosphere, the 

target planet's atmosphere and ionosphere, and the 
interplanetary plasma and solar corona. 

i. Atmosphere and Ionosphere of the Earth 

The neutral atmosphere iias an electrical path 
length of slightly under 3 meters in the zenith direction; 
this length increases as the secant of the zenith angle (co- 
secant of the elevation angle) . The lowest elevation 
angles usually employed in interplanetary radar experiments 
are about 6® for which the two-way effect on delay is 
about 0.2 ysec. Models for the. neutral atmosphere in current 
use are accurate to about 10% in their ability to predict the 
electrical path length at microwave frequencies. (This rather 
large uncertainty is due primarily to the difficulty in model- 
ling the very variable contribution of the atmospheric water 
vapor.) Thus, the uncertainty in the theoretical calculation 
of the contribution of the neutral atmosphere to the values of 
the time-delay measurements is never greater than 0.02 ysec. 
Moreover, for the Arecibo radar facility, the elevation angle 
can never drop below about 70° , implying a limit on our 
uncertainty in the atmosphere's contribution of under 0.002 ysec 
which will not be of any consequence. 

The effect of the earth's ionosphere is, of course, fre- 
quency dependent. For the radar frequency of 7840 MHz used 
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at the Haystack Observatory, the effect on the time-delay 
values has always been negligible compared to the measure- 
ment accuracy. At the 23S8 MHz radar frequency used at the 
Goldstone Tracking Station, the contribution of the ionosphere 
to the time delay during the day is of the same order as the 
contribution of the neutral atmosphere for high elev»ation 
angles; the ionospheric contribution, however, increases far 
less rapidly as the elevation drops below about 30°. Thus, 
the maximum contribution of the ionosphere to the delays 
measured at Goldstone is about 0.06 ysec and is small 
compared to their current delay measurement accuracy. 

At Arecibo, the currently used radar frequency is 
430 MHz for which the maximum contribution of the ionosphere 
to the delay measurements is about 0.6 ysec (recall the 70 
lower limit on elevation angles). This contribution can not 
at present be modelled with an uncertainty any smaller than 
about 30%, leaving the unmodelled part of the ionosphere con- 
tribution to the delay at about 0.2 psec for the daytime 
observations of Mercury and Venus. This uncertainty is 
almost as large as that accompanying the most accurate 
measurements currently being made in observations of these 
two planets at the Arecibo Observatory. For the past obser- 
vations of Mars near opposition, the signal passed through the 



-30- 


night-time ionosphere for which the ionospheric effect on 
delay is tenfold smaller, yielding a maximum uncertainty in 
the ionosphere contribution of about 0.02 usee which is neg- 
ligible. 

By the end of this year, before the next opposition of 
Mars, the upgrading of the Arecibo radar will be completed; 
subsequent planetary observations will be made at a radar 
frequency of 2380 MHz with a consequent 30-fold reduction in 

the influence of the ionosphere. 

i i . Atmo s phe re and Ionosph ere of the Target 
Planet 

Of the inner planets, only Venus has a 
significant neutral atmosphere from the point of view of af- 
fecting radar time-delay measurements. For observations of the 
delay to the subradar point*, only normal incidence is 
involved. Still, the incremental two-way delay introduced by 
Venus' neutral atmosphere is nearly 2 usee, virtually 


*The subradar point on a planet's surface is at the 
intersection there of the line drawn from the radar site to 
the planet's center of mass. 
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independent of the radar frequency*. Despite its relative 
enormity, this effect is not as devastating as it 
Sseems. Our only concern, in regard to relativity tests, 
is with the spatial or temporal changes in the atmospheric 
delay. The predominant contributions to the delay come 
from the lower few scale heights and, primarily because of the 
very large inertia of the lower atmosphere , temporal and 
spatial variations are expected to be no larger than a few 
percent. Indicative of the validity of this expectation is the 
finding from earth-based radio interferometry (Sinclair 
et al., 1972) that the variation of the 750°K surface 
temperature over the disc of Venus is under 15®K. Variations 
in the contribution of the atmosphere to the radar time- 
delay measurements that are caused directly by the topographic 
variations with respect to an equipotential surface can be 
"lumped" in with the topography effects (see Section III. 6) 
and need cause no additional concerns. It is possible, in 
principle, to separate the atmospheric from the strictly 


*The absorption of the radar signal in Venus' atmosphere 
is# however, a strong function of frequency, increasing 
approximately as the square. Thus the two-way attenuation 
for observations of the subradar point is about 7 db at X-band, but vir- 
tually negligible at the S-band frequencies used at Goldstone and soon 
to be used at Arecibo. 
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topographic effects on delay as discussed by Shapiro e t a l^ 

( 1973 ), but the implementation of the technique proposed in 
that article is not needed for relativity-test purposes. 

We note also that any error in modelling the 
(constant) delay introduced by Venus' atmosphere is of no 
concern; it is absorbed in the estimate from the radar data 
of the mean radius of Venus. Should an independent, and 
more accurate, estimate of the mean radius become available, 
as well as a better model for the atmospheric delay, then 
the comparison of tot radius, as nodified by the effect of the atmosphere, 
with the "radar" radius will allow another check to be made 
on the validity of the statistical filtering of the radar 

data . 

ionospheres on Venus and Mars have been measured using 
the radio-occultation technique (see, for example, Howard 
eti al. , 1974, and Kliore et al., 1972), but the integrated 
electron density in the zenith direction is, in each case, 
substantially smaller than that for the earth and poses 
no problem for the interpretation of either past or pro- 
posed radar observations of these planets. 

iii. Interplanetary Plasma and Solar Corona 

The sun emits, at high speed, 

charged particles that form the solar wind. This plasma, 
constituted primarily of electrons and protons, may for 
. purposes of discussion be divided into two regimes: the 
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interplanetary plasma and the "inner" solar corona; the 
demarcation line is arbitrary but we place it here at 
about four or five solar radii since in general only 
signals passing outside that limit can yield useful inter- 
planetary time-delay measurements. For signals that approach 
to about this distance from the sun, it becomes very dif- 
ficult to obtain accurate delay measurements at S-band radio 
frequencies because the turbulence of the corona severely 
restricts the coherence time of the radar signal. For such 
distances, even at X-band frequencies the turbulence ef- 
fects on the signals begin to affect the delay measurements 
adversely. At the 430 MHz radio fre- 

quency currently used by the Arecibo radar facility, useful 
measurements can not be made if the signal passes within 
about ten solar radii from the sun’s center. 

The above statements relate to average conditions of 
the solar wind. But the hallmark of the solar wind is 
its variability and, moreover, the unpredictability of 
this variability. It will be possible at times to make 
useful delay measurements at S-band for signals that approach 
within 4 radii of the sun; at other times, delay measurements 
will be thwarted for closest approach distances of 12 solar 
radii. Some examples of the enormous variability of the 
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turbulence of the solar wind are given in Part II of this 
report. Indications of the corresponding variability of 
the mean of the integrated electron density along paths 
passing near the sun are given in Figure 1 which is based 
on observations of the occultations of the Crab pulsar by 
the sun in 1973 (Weisberg et al . / 1975; see also Counse]itian and Rankin, 1972) 
Thus, for radar tests of relativity, one is restricted 
to signal paths that pass no closer to the sun than about 
four to five solar radii — unless new radar facilities 
are built which utilize substantially higher radio frequencies. 
But what about the interplanetary plasma beyond? How does 
it! affect the radar tests of relativity? Clearly the 
effect is greatest for the time-delay test (Section II. 2) 

Which depends on the signal path passing close to the sun. 

To obtain a quantitative estimate of the corruption of 
this test due to the interplanetary plasma we can utilize 
an approximate relation for the charged particle density, 

N(r), as a function of the distance, r, from the center of 
the sun: 

R 

N(r) ^ 5 X 10^ el/cm^; r ^ 5R , (3.1) 

r 0 

where R is the radius of the sun. This formula yields the 
I © 

average values for the charged-particle density near 

lOR but gives values somewhat too high near the earth's 
0 

orbit (r - 200 R ). The effect. Ax,, on an interplanetary 

0 P-*- 
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time-deiay measurement of such a plasma is given by (Shapiro, 
1964) : 


ATpi - (g^) [tan"^(^)+tan"^(^)l sec, (3.2) 

where f is the radar frequency in Hertz, d is the distance 

of closest approach of the signal to the sun's center, and 

X and X are the distances to the earth and target planet, 
e p 

respectively, from the point on the signal path closest 
to the sun. [When the earth and target planet are both on 
the same side of the sun, for example, near inferior con- 
junction or opposition, Equation (3.2) must be modified 
to indicate that the absolute value of the difference of 
the two arctangents must be evaluated.] For X^, Xp >> d, 
we obtain near superior conjunction 



0.9x10^^ 



R 

(-f)- 


sec , 


(3.3) 


showing that the plasma effect on time delay falls inversely 
with the impact parameter and with the square of the radar 
frequency. Since the predicted relativistic effect of 
solar gravity on signal delay, Ax^, given in Equation (2.2 ) , 
falls off logarithmically with d and is independent of 
frequency, the possibility for an accurate test clearly 
exists. For example, at the 7840 MHz frequency used 
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at Haystack/ =! 0.9 ysec for d = 5R^/ which is about 

0.5% of the corresponding predicted general relativistic 
effect of about 170 y sec. The signal-to-noise ratios 
achievable with the Haystack radar are, unfortunately, in- 
sufficient to take advantage of these favorable conditions. 

At S-band frequencies, doesn't even drop below 3% 

of At (with y = 1) until d =: 20 R for which distance 
r © 

- 3 ysec and Ax^ - 115 ysec. This "cross-over" point 
is illustrated in Figure 2 for the S-band radar frequency 
of 2380 MHz to be used by the upgraded Arecibo facility 
with Venus the target planet. 

Because of the use of S-band at Arecibo, nowhere 
near the full advantage of the available signal-to-noise 
ratio can be taken in the performance of the time-delay, 
or signal-retardation, test of general relativity. For 
the other relativity tests which depend only on the orbital 
motion of the inner planets the situation is much ameliorated 
but the effects are still not negligible. For example, at 
the elongation of Mercury, Equation (3.2) shows that for 
f = 2380 MHz, ysec, whereas 

at inferior conjunction 0*1 usec. For 

Venus, the corresponding elongation and inferior conjunction 
numbers are slightly under 0.15 ysec and 0.04 ysec, res- 
pectively. For Mars at opposition, the most favorable case, 

we obtain Ax , 0.03 ysec. The above results are summarized 

pi 

in Table 1. 
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Table 1. Average time delay, in ysec, expected from inter- 
planetary plasma for 2380-MHz radar observations 
at different orbital configurations. 


Planet 

Mercury 

Venus 


Elongation 

0.3 

0.15 


Inf. Conjunction 
0.1 
0.04 


Opposition 


Mars 


0.03 
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These estimates of the effects on delay of the 
interplanetary plasma between elongations and inferior con- 
junctions or oppositions are apt to be , on average , nearly 
a factor of two too large since Equation (3.1) yields a 

3 

value for N of nearly 11 el/cm at 1 a.u., whereas the 

3 

average value is perhaps only about 5 to 7 el/cm . Fur- 
thermore, using the plasma data obtained from the Pioneers 
and other relevant instruments, we should be able to model 
the plasma effects on delay with an uncertainty perhaps as 
low as 50%. 

The above discussion demonstrates that the analysis and 
initerpretation of radar data obtained in the past, with the 
exception of some Arecibo observations , have not been serious 
ly compromised by the interplanetary plasma. Arecibo and, 
perhaps, Goldstone observations to be made in the future 
may be seriously affected. The most elegant way to virtually 
eliminate this corrupting effect is to make radar observa- 
tions simultaneously, or nearly simultaneously, at two 
frequencies. At Goldstone there is now the capability 
to make observations at both S-band and at X-band. At the 
much more powerful Arecibo radar facility, only an S-band 
capability is currently planned. However, should the re- 
surfaced antenna at Arecibo be sufficiently well adjusted 
to be usable efficiently at X-band frequencies, an X-band 
radar system could be installed there too* 
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Delays measured simultaneously at two radar frequencies 
can be combined to yield a delay observable freed from 
•plasma effects. The standard deviation, a(x^), of this 
combined observable is related to the standard deviations 
a(T^), i = 1, 2, of the delays measured at the individual 
frequencies f^ by: 


cjCt^) {■ 


fy (^ 2 ^ 


2,2 


(f -f ) 
'^2 1 ^ 




1/2 


(3.4) 


For fj^^2 GHz and £2 - 8 GHz, Equation (3.4) yields 


a(x^) 1.07 cf(T^), (3.5) 

for a(Tj^) = a(T 2 ). Thus the "penalty" paid in freeing 
the delay measurements from' plasma effects is only about 
7% for a dual, S- and X-band system. 

For the purposes of the relativity tests, with X-band 
available, the S-band capability is not required except 
near superior conjunction where its effectiveness will 
cease only when the turbulence in the corona causes the 
delay measurement accuracy to decrease significantly. From 
elongation through inferior conjx’.nction, even for Mercury, 
the maximum effect of the interplanetary plasma at X-band 
will be - 0*02 ysec which is nearly negligible. 
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4. Gravitational Perturbations 

The gravitational perturbations to the target 
planet and to the observing platform come from many 
sources: the galaxy, neighboring stars, the sun, the planets, 

asteroids, comets, and satellites. We treat each of these 
briefly, in the order mentioned, 
i. Galaxy 

The galaxy introduces a gradient in the gravi- 
tational field so that the bodies in the solar system will 
thereby experience differential accelerations. For the pur- 
poses of an order-of-magnitude estimate we may consider the 
entire mass of the galaxy to be concentrated at its center. 

The maonitude of the differential acceleration, a , of a planet relative to 
^ • y 

the sun caused by the galactic mass will thus be given ap- 
proximately by 



(3.6) 


where G is the constant of gravitation, the mass of the 
galaxy, r^ the distance from the galactic center to the 
sun, and A the distance between the planet and the sun. 

Using the approximate values (see, for example, Allen, 1963) 
M =: 10^^ M = 2 X 10^^ gm and r 9 kpc 3 x 1(?^ cm, we 

g © y 

obtain for Mars (A == 2 x 10^^ cm) , the outermost planet of 
interest, the value 

Sig - 2 x lo”^^ cm/sec^ (3.7) 

since G ^ 6.1 k 10~® in cgs units. Since the most im- 
portant, "in-plane” effect of this perturbing acceleration 
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will tend to average out over an orbital period, we can 

obtain a crude approximation to the magnitude of the effect 

2 . 

by using the familiar formula s =(1/2) at with t being 
one-quarter {- 1.5 x 10^ sec) of Mars' orbital period; we 

obtain 

s =: 2 X 1 0 ^ cm , (3.8) 

g 

which is utterly negligible. The gradient in the gravi- 
tational field due to neighboring galaxies is obviously far, 
far smaller. 

ii. Neighboring Stars 

Apart from the overall effect of the 

galaxy, one should, for completeness, estimate separately 

the influence of stars in the solar neighborhood. Using 

Equation (3.6) with Mg replaced by M^ and rg by r^, and with 

1 fl 

M M and r == 1 pc 3 x 10 cm, we obtain for Mars: 
s © s 

a - 2 X lo”^^ cm/sec^; s =: 2 x 10 ^ cm, (3.9) 

S s 

again completely negligible. Thus, to this extraordinary 
degree, the solar system can be treated as a closed 
dynamical system. 

iii. Sun 

Here we are concerned with the gravi- 
tational effects of the sun on the motions of the planets, 
apart from the direct Newtonian and relativistic influences 
of an equivalent point mass. Three such effects 

are of potential significance: the mass distx'ibution within 
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the suii/ changes in the mass of the sun/ and tidal bulges 
in the sun. We discuss each in turn. 

iii. 1. Mass Distribution 

The nonspherically symmetric 

part of the mass distribution of the sun has been the sub- 
ject of considerable controversy over the past seven years. 

Dicke and Goldenberg (1967) claim to have observed a 
visual oblateness of the sun and interpreted it, using 
Von Zeipel's theorem (see, for exaTiple, Clayton, 1963), to jjiply that 
the coefficient of the second harmonic of the sun's gravitational field 

had a value about two orders of magnitude larger than would 
be expected were the sun to be rotating uniformly with depth 
at the angular velocity observed at the surface. In par- 
ticular, Dicke (1974) concludes that 


J_ (2. 4 ± 0. 3) X 10 
2 © 


(3.10) 


where Jo is the coefficient of the second-degree zonal harmonic in the usual 

. 2 Q) 

spherical harmonic representation of the gravitational potential: 

GM ^(5)2 i 

V(r,e,<(.) = - ^ {1-^2 


where R is the equatorial radius of the sun; r, 6, and 
0 

(j) are the radius, colatitude, and longitude of the 
(external) field point; and P 2 is the second-degree 

Legendre polynomial. [See, for exartple, Dwight (1949) for the normalization 
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used . ] We have estimated the value for J 2 ^ expected on 
the basis of uniform rotation and find 


J . = 1 X 10 

2 © 


(3.12) 


with the uncertainty, even given uniform rotation, being 
perhaps as large as a factor of two. 

Dicke ' s interpretation of the Princeton solar- 
oblateness measurements have been criticized by many. 
With the proper interpretation of the Princeton data 
in doubt, it is clear that other experiments are 
required. Older data, obtained with a heliometer over 
y 0 ar p©riod by Schur and Ambronn at the turn 
of the century, indicated the absence of any visual 
oblateness of the sun with an uncertainty several- 


fold below the magnitude of the effect reported by Dicke 
and Goldenberg. But Dicke has raised the possibility that 
systematic errors accompanying the use of the heliometers 
may have significantly affected this apparently null 
result. More recently, however. Hill et al. 
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(1974) using a technique very different from that employed 
at .Princeton to define the solar limb, have also reported essentially 
a null result, supporting the earlier conclusions of Schur 

and Ambronn. 

Rggardiess of the final outcome of the controversy 
surrounding the direct measurements of the visual oblate- 
ness of the sun, it is only the gravitational oblateness 
which is of. concern for the radar tests of relativity. 

Since the gravitational oblateness characterized by J 2 ^ 
can in any event only be inferred from the visual oblate- 
ness in a theory-dependent (and somewhat controversial) man- 
ner, we would like to be able to determine J 2 ^ directly 
from its dynamical effects. The perturbations of the 
planetary motions due to J 2 ^ are unfortunately smaller than 
those due to y and 3 (see Section II) . Even with J 2 ^ as 
large as 2.5 x lO”^, the perturbations of Mercury's orbit are 
only about 10% of those due to y and 3. With interplanetary 
time-delay measurements, the effects of J 2 ^ can be separated 
from the relativistic ones due to the differences in 
(1) the dependence on orbital radius of the secular advance 
of the perihelion (the J 2 Q effect falls off faster by one 
power of the inverse radius ) ; (2) the short-period pertur- 

bations of the orbits; and (3) the relations between or- 
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bital radius and period for a suite of planets. The rela- 
tive efficacy of these differences in effecting a separation 
has never been studied with care. In particular, the 
‘relative importance for these differences of the "masking" 
effects of the orbital elements, which must be estimated simul- 
taneously with >^20' never been evaluated. But 

even after the relative efficacies have been determined, and 
a suitable strategy of observations adopted, and even if J 2 ^ 
is of the order that one would predict from uniform rotation, 
the necessity for the precise determination of J 2 ^ from its 
dynamical effects will likely remain the nost iirportant limitation 
so far mentioned on the accuracy achievable in the radar 
tests of relativity that involve planetary motions directly. 

The third and higher-order harmonics of the sun's 
gravitational potential are obviously of no concern for the 
study of planetary motions. In addition to the overwhelming 
likelihood, based on theoretical considerations, that the 
coefficients will be much smaller than , their effects fall 
off with increasingly higher powers of the ratio of the 
sun's radius to the planet's orbital radius. Thus, even 


were of the same order as the effect on Mercury 

30 -^0 

orbit would be almost two orders of magnitude smaller. 



iii. 2. Mass Changes 


The second effect of the sun, 
mentioned at the start of this subsection, concerns its 
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change of mass. The main mechanisms for mass loss are 
electromagnetic radiation and the solar wind; the effects 
of neutrino radiation are negligible by comparison. Mass 
accretion takes place through direct impact of (perturbed) 
orbiting bodies and through the decay of the orbits of dust 
particles under the influence of the Poynting-Robertson ef- 
fect (see, for example, Robertson, 1937) . The losses can 
be evaluated more accurately than the gains. For electro- 
magnetic radiation, we have 


dM 

O 

dt 


radiation 


-4TTr^a, — -4x10^^ gm/sec 

-6 X lO"^'* MVyr, 

© 


(3.13) 


where r 1.5 x 10^^ cm is the radius of the earth's 
orbit; - 2.0 cal/cm^-min = 1.4 x 10^ ergs/cm -sec is the 

flux from the sun at 1 a.u.; and c 3 x 10^^ cm/sec is the 
speed of light. We assume, of course, that the radiation 
is spherically symmetric. We make the same assumption, 
with substantially less reliability, for the particle flux 
in order to get a global estimate: 



=: -4‘TTr^nmv - -2 x 10 gm/sec 

particles 

= -3 X 10“^^M^/yr , 



(3.14) 
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where n = 10 cm”^ is the (over) estimated average proton number 
density at 1 a.u.; m - 1.67 x 10 gm is the mass of a 
proton; and v ^ 400 km/sec is the (over) estimated average 
velocity of the solar wind at 1 a.u. Particles other than 
protons make a much smaller contribution to the mass loss; 


their omission from dM /dt is far less serious than the 

© 

uncertainty in the estimate of the proton mass loss. 

The mass accretion is difficult to estimate reliably. 


Conceivably, but very unlikely, it might even exceed 

lO”^^ M yr”^ and lead to a net mass gain for the sun. 
© 

However current estimates (Marsden, 1974) 


indicate a value for dM /dt 

0 


. of about 10 M yr ^ 

influx © 


or smaller. For purposes of discussion let us assume that 

the change in mass of the sun, whatever its sign, doesn't 

exceed lo”^^ M yr”^. Such a bound implies that the uncer- 
0 

tainty in the mass changes of the sun will not influence 


any of the radar tests of relativity for many years to come. 
The most pronounced effect of these mass changes will be a 
slow, outward spiralling of the planetary orbits which would 
affect most importantly the detection of any possible change 
in the gravitational constant. But even for radar time- 
delay measurements spanning two decades, the sun’s mass 
changes, if within the limits specified, would introduce 
a maximum change in a delay measurement of under 0.4 ysec. 
For the near term, therefore, solar mass changes will not 
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limit the accuracy to which one can discern any possible 
change in the gravitational constant. 

iii. 3. Tides 

Tides are raised in the sun by 
the planets. Such tidal bulges affect the planetary orbits / 
in particular the semimajor axis, the eccentricity, and the 
inclination*. However, the magnitude of these effects is 
negligible because of its dependence on both the square of 
the (small) planetary mass and the inverse sixth power of 
the ratio of the planet's distance to the sun's radius. This 
tidal effect on planetary orbits is also inversely propor- 
tional to the relevant "Q" of the sun. But even for Q - 10 , 
likely a gross underestimate, and for Mercury, the most 
strongly affected planet, we find that the major effect, on 
Mercury's semimajor axis, is smaller than that due to the 
sun's mass loss and so may be neglected for the foreseeable 
future in regard to radar tests of relativity. 

iv. Planets 

The planets affect the radar time- 
delay measurements because of the perturbations they cause 
in the orbits of the earth and target planet. These effects, 
of course, depend on the masses and orbits of the perturb- 
ing planets. We discuss the orbits first. A third effect, 

*Because the planets are all beyond the radius of synchronous 
rotation with the sun, the orbits will tend to spiral out; they 
will also tend to circularize. 
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due to the spin-orbit coupling, will be discussed last. 

iv. 1. Orbits 

For the inner planets , where direct 
radar measurements exist for each, it is necessary to es- 
timate simultaneously with the relevant parameters for the 
relativity tests the six orbital elements for each inner 
planet . These estimates must be based on the available radar and 
optical data concerning these bodies. It is desirable, but not mandatoiry , to 
include simultaneously in the analysis the optical data for the outer 
planets and the parameters describing their orbits (and masses) . 
This expansion of the data and parameter sets in the solution 
is not now mandatory because the outer planets are of interest 
•to the relativity tests under discussion only insofar as 
they perturb the orbits of the inner planets. But an error 
in an element of the orbit of an outer planet has only a 
second-order effect on the determination of the orbit of an 
inner planet. Nonetheless, the fractional uncertainties 
accompanying the radar time-delay observations of the inner 
planets are becoming so much smaller than those for the op- 
tical observations of the outer planets that, for example, 
the initial longitude of Jupiter in its orbit will soon be 
able to be estimated with higher accuracy from radar ob- 
servations of Mars than from the existing optical observations of 
Jupiter. At this point the weak coupling between the inner 
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• and outer planets will become sufficiently important to 
require a simultaneous solution in order to obtain an 
accuracy for the relativity tests commensurate with the 
inherent accuracy of the radar data. 

Before this "crossover" point is reached, however, it may be possible, 
with the upgraded Arecibo radar facility to 

observe the Galilean satellites of Jupiter directly and 
possibly, also, Saturn's largest satellite. Titan. From 
the analysis of these direct radar measurements of the 
satellites of Jupiter and Saturn the orbits of the latter 
two will be determined far more precisely than would be 
possible from the inner-planet radar data. On the other 
hand, the inner-planet radar time-delay data, because of the 
much shorter distances involved and the inverse fourth 
power dependence of the echo strength on. distance, will remain of 
far higher accuracy than the corresponding data from the 
outer solar system and so the latter data would add little 
to the determination of the orbit of the earth — the cannon 

observing platform. This last statement also leads to the 
cponclusion that, with no appreciable loss of accuracy, the 
inner-planet data can be analyzed separately from those of 
the outer planets for radar tests of relativity. 


iv. 2. Masses 

Several comments are required in 
regard to planetary masses. For the inner planets, radio- 
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tracking data from spacecraft near planetary encounters 
have yielded the masses for Venus, the earth, and Mars with 
‘uncertainties at the level of several parts per million 
(for a summary, see Ash, Shapiro, and Smith, 1971), Pre- 
liminary analysis of the Mariner 10 data (Howard et al, , 

1974) give an uncertainty for Mercury's mass of about 

4 

1 part in 10 . The radar time-delay data also contain in- 
formation on these masses due to the perturbations induced 
by these planets in one another's orbits. Since these planets do 
not make close approaches, the spacecraft data are intrin- 
sically more powerful for the estimation of planetary masses 
and all but eliminate the need to estimate them from the 
radar data. However, inclusion of the mass parameters in 
the analysis of the radar data for the relativity tests 
provides an important check on the existence of unsuspected 
S'- sterna tic errors. Any disagreements between the "radar" 
and "spacecraft" values for the masses, large compared to the 
accompanying standard deviations, would imply the presence 
of such errors. 

The outer planet masses are best determined from mutual 
perturbations among these planets, with two exceptions : 

Jupiter and Pluto. The mass of Jupiter (plus its satellites) 
is at present determined most accurately from analysis of 
asteroids whose mean motions are commensurable with Jupiter's. 
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The Pioneer 10 flyby of Jupiter provides the next most 
accurate value. For Pluto, despite some claims to the 
contrary (see, for example, Seidelmann et al_._, 1971), we 
have shown that its mass can not be determined reliably 
from its perturbations of the orbits of the other outer 
planets, Pluto’s mass is simply too small and the optical 
observations of the outer planets are too crude to allow a 
determination (see Ash, Shapiro, and Smith, 1971). By the 
same to).en, Pluto is not massive enough to affect the radar 

tests of relativity. 

AS for the inner planets, however, the estimation of the 
masses of the outer planets from the inner-planet time-delay 
data provides a useful cheg’t. And, in fact, the formal ac- 
curacy achieved in such estimates of the masses of Jupiter and 
Saturn are only several-fold poorer than for the best now avail- 

able. 

iv. 3. Spin-Orbit Coupling, 

The spin-orbit resonance of Mercury 
(see, for example, Colombo and Shapiro, 1966, Goldreioh and 
Peale, 1966, Counselman, 1968) is potentially of direct signi 
ficance for radar tests of relativity. This resonance is 
responsible for a slight additional advance, of Mercury’s 


perihelion given by 


, 21TT .B-A. 26 


rad/rev. 


(3.15) 
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where C are Mercury's principal moments of 

inertia; a, e, and are Mercury 's orbital semi-major axis, ec- 
centricity, and radius, respectively; and 6^ is the angle betv^een Mer- 
cury's axis of miniiram moment of inertia and orbit major axis at perihelion. 
Since the value of (B-A)/C is unknown for Mercury, the pre- 
cise magnitude of this effect can not be determined. 

However, based on analogy with the moon and on the likely 
extent to which Mercury could remain out of hydrostatic 
equilibrium over the relevant time interval, we conclude 
that (B-A) /G <2 X lo”^. This bound on the fractional 
difference in the principal moments of inertia irtplies a peri- 
helion advance of no greater than 0'.'003 per century. Hovever, in order to 

improve the corresponding test of general relativity to the 

4 

point where the uncertainty is a part in 10 or less/ it 
will be important to have an independent estimate of (B-A)/C 
of uncertainty under 2 x lo”'^. Further analysis of the 
Mariner 10 radio tracking data obtained near Mercury encounters 
may yield a determination of sufficient accuracy; the cor- 
responding data from an orbiter certainly will. 

V. • Asteroids 

The asteroids, or minor planets, 
pose the same generic problems as the major planets. But 
because of their large number, their location in the solar 
system, and their small masses, the detailed problems posed 
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are quite distinct. A considerable effort has been expended 
under our direction (see Friedman, 1970, and Sherman, 1973) 

’to’ develop the most feasible approach to include the aster- 
oids in the data analysis. One cannot simply follow the 
"ideal" procedure and include all the asteroids in a grand 
and glorious N-body integration of the motion of the solar 
system bodies, with the associated orbital elements and masses 
added to the parameter set: N is just too large and, in 
fact, is not even known (except as a lower bound) . As a 
reasonable compromise with practicality, we have chosen the 
apparently twelve most massive asteroids (see Allen, 1963, 
and also Matson, 1971, 

regarding the proper place of Bamberga) , and we have taken 
their orbits as given by the separate reductions of the 
optical observations of these objects. Thus, in our in- 
tegration of the planetary orbits, we include the effects 
of these asteroids on the planets, but ignore any possible 
corrections to the asteroid orbits that are implied by the 
(new) planetary orbits. This procedure is of more than suf- 
ficient accuracy for our purposes; the same arguments ad- 
duced for the major planets apply here, but a fortiori , since 
the proximity of the asteroids to the inner planets is, 
except possibly for a "pathological" case of near encounter, 
more than offset by the far smaller masses of the asteroids. 
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The masses of these dozen asteroids are set in accord with 
other determinations, based mostly on asteroid-asteroid 
perturbations, or estimates, based on measured diameters ^d 

3 

an assumed density of 3 gm/cm . We also have the capa- 
bility to add any subset of these masses to the list of 
parameters to be estimated in a given weighted— least- 
squares analysis of the data. The current limit of twelve 
could easily be raised to the order of thirty were such a 
change warranted. 

What of the myriad of other asteroids? How can their 
effects be incorporated? We have chosen to replace the 
remainder of the asteroids by a circular mass ring whose 
radius, mass, inclination to the ecliptic, and longitude of the node 
can be either fixed or estimated as desired in any given analysis. 

The effects of the statistical fluctuations of the true 
asteroid distribution about this average ring must be 
evaluated in order to ascertain the usefulness of the ring 
model. {For a preliminary analysis, see Sherman-, 19 73.) 

Finally, we comment on the magnitude of the asteroidal 
perturbations on the orbits of the inner planets. The 
largest asteroid, Ceres, has a mass of approximately 
6 X 10”^° M (Schubart, 1974), which may represent about 
a third* of the total mass of the asteroid belt. We can 


*This is the fraction inferred from Allen (1963) , but it 
may be very unreliable. 
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make a crude estimate of the maximum effect its perturbations 
will have on inner-planet time-delay data, by evaluating 

• the two-way light-time equivalent, ^'^^eres' dist 

that Mars would be perturbed during a quarter of its or- 
bital period if the acceleration due to Ceres acted in a 
constant direction with a magnitude corresponding to the 
distance of closest approach of the two objects (see Sub- 
section III.4.,i). We find: 0.7 ysec. A more 

accurate, but still imprecise, bound obtained by comparison 
of two numerical integrations of Mars' orbit, one including 
and one omitting the effect of ' Ceres , yields a compatible 
bound.* These numerical values indicate that the fractional 
error in the estimate of an asteroidal mass need be no less than 

• about 0.1, except for the case of Ceres, in order to insure 
that, with the expected capability of the Arecibo radar 
system, no significant corruption of the radar tests of 
relativity will result, it is also clear that, for the same pur- 
poses, the gravitational effects of ccmets can be ignored. The corres- 
ponding effects of interplanetary "dust” have not been estimated expli 

citly but are presumably negligible. 

vi. Satellites 

The corrupting effects of satel- 
lites are nost easily discussed by separating them into three classes : 
the iToon, Phobos and Deimos, and the satellites of the outer planets. 


‘This technique f "?=^erencr!n°tota? masf in\he 

♦■he r>iar>pt's lonqitude due to the dirrerence xn xvjv-a 

two ^aSuLLons? this drift would largely disappear in an 
"operational" comparison. 
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vi . 1 . The Moon 

The moon affects interplanetary 
time-delay measurements primarily by causing the earth to 
revolve about the earth -moon barycenter. However, any in- 
accuracies in our knowledge of the geocentric orbit of the 
moon are reduced by a factor of 81 — the earth-moon mass 
ratio — in their effect on our calculation of the bary- 
centric motion of the earth. Because of this nearly two 
order of magnitude reduction factor, only the lack of suf- 
ficient precision in knowledge of the earth-moon mass ratio 
and, to a lesser extent, in the moon's initial mean anomaly 
has been of any concern in the analysis of the inter- 

planetary radar data. Thus these two parameters are in- 
cluded in the set to be estimated. 

In the future, laser observations to the retroref lectors 
oh the moon -- currently being made with errors of the order 
of only a few nanoseconds — will yield an orbit that 
more than ireets the accuracy requirements for the calculation 
of the earth's barycentric motion in the 

analysis of the radar data. The perturbations induced by 
the moon in the orbits of the other inner planets can, of 
course, be calculated with far higher accuracy than will 
be required by the interplanetary time-delay data. 

vi . 2 . Phobos and Deimos 

Aside from our moon, Phobos and 

Deimos are the only known satellites of the inner planets. 
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These two satellites of Mars will, of course, cause Mars to 
rotate about the corresponding barycenter. The orbital 
radius of Phobos is only about 40% of that of Deimos, but 
the former's approximately five times larger volume (Dux- 
bury, 1974) makes it, under the assumption of equal 
density, the greater contributor to the 

amplitude of Mars' barycentric motion. Assuming the density 
to be 3 gra/cm^, we find that the amplitude of Mars' barycentric 

motion will be 
M 

At, - - a - 0.002 usee (3.16) 

b c M^ P 

where M |-tt abc p 1.8x10^^ gm is Phobos' mass; 
p o 

33 3 

V 3b9§7 W ^ m^ssi ap 9x10 km 

is the semimajor axis of Phobos orbit (Allen, 1963) ; and 

a - 13.5 km, b 11.5 km, c - 9 . 5 km are the principal 

axes of Phobos (Duxbury, 1974), with p its mean density. 

Thus, for radar tests of relativity. Mars' barycentric motion 

is irrelevant. 

We might also remark that any as yet undetected satel- 
lites of either Mercury or Venus would be too small to be 
relevant. For the likely limit on their size, such 
satellites could not be far enough from their parent planet 
to produce a noticeable effect, because orbits of that 
size would be unstable against solar perturbations. 
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vi . 3 . Satellites of Outer Planets 

The outer planets affect 

,th.e presently contemplated radar tests of relativity only 
by their perturbations on the inner planets. To bound the 
corresponding effects of the satellites of the outer 
planetS/ we need only consider Callisto's effect on Mars* 
orbit. Because Callisto is in close orbit about Jupiter, 
the former's effect on Mars would be essentially quadrupolar . 
We can approximate the instantaneous acceleration exerted 
by Callisto on Mars by the simple, easily derived expression: 

a 3GM ■ lO"^^ cm/sec^/ (3.17) 

c c ^4 

^J-M 

11 

where M 10'^° gm is Callisto's mass and ~ 2 x 10 cm 

C t 

its mean distance from Jupiter (Allen, 1963) , and where 

r, 3.7 a.u. - 5.5 x 10^^ cm is the distance of closest 

J-M 

approach of Jupiter and Mars. Considering Mars' orbital 
period, we find that the maximum effect on an earth- 
Mars time-delay measurement would be: 


c c c 4 


(■ 


-4 

-7x10 usee. 


(3.18) 


where denotes Mars' orbital period. Again we conclude 
that the problem caused by satellites is negligible. But 
in any event, we could calculate their effects to high accuracy so 

there is really no problem. 

5. Non-Gravitaticnal Perturbations 

we consider now the direct effects of 
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nongravitational accelerations on the orbits of the inner 
planets. These come primarily from three sources: electro- 

magnetic radiation, the solar wind, and particulate matter. 

The third category covers impacts on the planet of inter- 
planetary dust, cometary material, and asteroids. We es- 
timate the effects of each in this section. 

i . Electromagnetic Radiation 

The acceleration of a planet due to 
electromagnetic radiation {sunlight pressure) 
is proportional to the area-to-mass ratio of the planet and 
decreases inversely with the square of the distance to the 
sun. Of the planets of direct interest to us. Mercury is the 
nost severely influenced by this perturbation and will be used to place 
an upper bound on it. Vfe first consider the radial conponent of the pres- 
sure which produces an acceleration on Mercury equivalent 
to a slight decrease , AM , in the mass of the sun. 


We may bound this equivalent decrease by 


(^) 4 

where (A/M)y 5.6 x 10 
ratio for Mercury*; (I^/c) 


•5 1 . 8 X 10^® gm, (3.19) 

2 

cm /gm is the area-to-mass 
-5 2 

- 4.65 X 10 dynes/cm is the 


By area, of course, we mean the projected area. 
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13 

pressure of sunlight at 1 a.u., and ~ 1.5 x 10 cm is 
1 a.u., the semimajor axis of the earth's orbit. The 
factor of two on the right side of Equation C3.19) re- 
flects the overestimate (upper bound) based on conditions 
pf "flat-plate" reflection. The magnitude of the effect on 
time-delay measurements from earth to Mercury of this 
effective decrease in the solar mass, a radial acceleration, 
can be bounded by 


At I 5 X 10 ^ tysec, 

sr c b M 0 Py 

>13 


(3.20) 


where r^ = 0.6 x 10^^ cm is the semimajor axis of Mercury's 

¥ 

orbit; - 0.24 years is Mercury's orbital period; and 
t is expressed in years. (Note that the factor of two in-, 
troduced by round-trip delay is cancelled by a factor of 
one half introduced by the relation between the 


sensitivity of the rate of change of longitude displacement 
to the effc^ctive change in solar mass.) For radar observa- 
tions spanning several decades, this sunlight-pressure 
perturbation can be ignored . 

There is the possibility that the reflection properties 
of Ifercury, for example, could be sufficiently "skewed" for the 
reflected radiation to cause the average acceleration vector to 
have a substantial fractional component in the tangential 
direction that is parallel to Mercury's direction of motion. 

(All surface longitudes on ^tercury do not receive equal amounts of solar il- 
lumination, either in magnitude or direction, because of the spin-orbit 
resonance and the orbital eccentricity. Ihus, even if the surface reflected 
isotropically, longitude variations of albedo could result in a nonzero- 
average tangential radiation. ) 
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Displacements due to such accelerations grow, of course, 
qviadratically with time. Thus a bound similar to that 
given by Equation (3.20) can easily be derived, for the tan- 
gential acceleration: 

i ^ t^ 1.1 X lO"^ F t^ ysec , (3.21) 

' 

where F is the fraction representing the ratio of the mag- 
nitude of the tangential acceleration due to sunlight pres- 
sure to that of the total due to sunlight pressure. The 
time, t, is again given in years. For F of the order of 
0.01, which is much higher than crude estimates would indi- 
cate, a set of radar observations spanning more than a decade 
might begin to be noticeably affected by such a tangential 
acceleration. 

A related source of a tangential acceleration is the 
infrared radiation from a planet with a temperature dis- 
tribution asymmetric with respect to the sun-planet line. 

We can place a crude upper bound on this acceleration with 
the following easily-derived formula: 


a . 

ir 


3aT AT ,A, ^ , , ..-14 . 2 

— c 3.4 X 10 cm/sec 


(3.22) 


^ c O A 

where a 5.67 x 10 ergs/cm -deg K -sec is the Stefan- 
Boltzmann constant, T 600°K is an upper limit on the 
average surface temperature on Mercury, AT 50®K is an 
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upper limit on the average difference in temperature between 
the "leading" and "trailing" hemispheres of Mercury; and, 
as stated above, (A/M)^ = 5.6 x 10 cm /gm is the area- 
to-mass ratio for Mercury. In this approximation, Mercury 
is being treated somewhat like two flat plates whose normals 
are in the direction of orbital motion and whose areas are 
equal to the projected area of Mercury. 

The ■ maximum effect on the earth— Mercury time— delays of 
this "infrared" acceleration is given by: 

At. - a. t^ =; O.OOlt^ ysec , (3.23) 

ir c ir 

where t is measured in years. If this upper bound actually 
j- 0 pj- 0 sented the effect, we can see that the interpretation 
of a span of data more than a decade in extent would begin 
to be seriously affected by this infrared radiation. Thus, 
a more careful bound on this effect should be developed 
before too long a span of data is analyzed. 

Finally, we note that the Poynting-Robertson effect 
(Robertson, 1937) , mentioned earlier, also introduces a 
tangential acceleration but smaller in magnitude by a factor 
of about v/c relative to the radial acceleration due to 


sunlight pressure. Thus, 
V 

a 


“ 5 X cm/sec^ 

pre c M 9 c 




(3.24) 


where v^ is the orbital velocity of Mercury and ir^/r^) is 
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the ratio of the earth's orbital semi-major axis to Mer- 
cury's. This acceleration is clearly negligible compared 
to the upper bound calculated above for t±ie infrared radiation. The or- 
bital effects of the cosmic background black-body radiation are also negligible. 


ii . Solar Wind 

To obtain approximate 


values for the effects of the solar wind on planetary orbits, 
we follow the same procedures as used in Subsection 4. iii.2. 
If we assume that the solar wind impinging on a planet's 
surface is completely absorbed (we can ignore the effects of 


the resultant increase in mass of the planet) , then we can 


approximate the radial acceleration due to the wind by 




-17 , 2 

6 X 10 cm/sec , 


(3.25) 


where (nmv^)^ represents the radial momentum transfer to 
the planet per unit area per unit time; and where, indi- 


vidually, n ~ 70’ particles/cm^ is the density of protons 
in, and v - 300 km/sec the speed of, the solar wind at 
Mercury's orbit with m = 1.67 x gm being the mass of a 

proton. To the extent that (nmv^) might be constant over 
Mercury's orbit, the effect of a^^. on the orbit would be of 
short period only. In this case, we can place an upper bound 
on the consequent effect on time delays by using the 
"s = ^ at^" method with t one quarter of Mercury's orbital 
period. We find Ax 10”%sec, which is completely neg- 
ligible. To the extent that nmv^ varies inversely with the 
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square of Mercury's distance from the sun we have a bound 
similar to that given by Equation (5.2), except that the 
wind's effect is smaller by a factor 


(nmv^) 


wr 




2 X 10“^, 


2 (^) (^) 2 

? . K 

where values for the various quantities were given above. 


(3.26) 


For completeness we also point out that, due to aber- 
ration , the solar wind causes a tangential acceleration 
as well, given approximately by: 

a ^ ^ -Si a =il.5x lO”^^ cm/sec^ / (3.27) 

wt v„ wr 
w 

where v^ = 50 km/sec is Mercury's orbital velocity and 

v 300 km/sec is the velocity of the solar wind at Mercury's 

w 

orbit. Comparison with Equation (3.24) shows that a^^ is about 
a factor of three smaller than the corresponding acceleration 
due tOi the (negligible) Poyntinq-Robertson effect. 

Althou^ no explicit bounds have been calculated, it is probably 
quite safe to conclude that the effects of the magnetic field and of the 
particles in the solar wind on the planetary magnetospheres (where they 
exist) introduce no discernible effects on the orbits of any of the inner 
planets. Similarly, we conclude that the effects of cosmic rays are negligible. 

iii. Particulate Matter 

It is difficult to estimate 

reliably the momentum transfer and mass changes that accrue 
from collisions between planets and interplanetary dust. 


- 68 - 


comets, and asteroids. According to Marsden (1974), 

the influx of such material on the earth at present probably 

12 

averages no more than about 5 x 10 gm/yr. The consequent 
fractional mass change in the earth of 10 yr ^ will have 

negligible import for tests of relativity as can be seen 
from analogy with the discussion in Subsection 4.iii*2. 
momentum transfer to the earth is undoubtedly a minuscule 
fraction of the total of the magnitudes of the momenta of 
.the individual particles. However, even if the net transfer 
were equal to this total, the effect on the earth's orbit 
would be small. Taking 20 km/sec as the maximum relative 
velocity, v^., and assuming all these velocities to be parallel 
to the earth's orbital velocity, we obtain a gross upper bound 


On the tangential acceleration due to particulate matter of 

(3.20) 




“17 2 

5 X 10 cm/sec , 


where = 1.7 x 10^ gm/sec is an upper bound on the rate 
of mass gain by the earth as estimated above. Comparison 
with the previous subsection shows that even if this upper 


bound were the actual tangential acceleration, the latter 
would have negligible effects on the relativity tests. We 
may also reasonably conclude that the influx of particulate 
matter on the other inner planets can also be neglected. 

Only Mars is likely to have a substantially higher influx 
than the earth, but not by a factor large enough to render 
invalid the above conclusion. 

Any, perforce infrequent, collisions between inner planets 

18 

and large objects (~ 10 km in diameter and 1.5 x 10 g in 
mass) could be modelled as impulses in the analysis. 
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6 . Planetary Topography 

P^^haps the most vexatious problem hindering 
tile maximum utilization of radar dats for testing general 
relativity is that caused by the large topographic 
<rariations over the surfaces of the target planets. Mars, 
the worst offender, introduces variations in the time 
delays of up to about ± 50 ysec. Fortunately, these topo- 
graphic variations over the surface are constant in time, 
at least on scales of current interest for radar tests of 
relativity. Still, we can not simply conclude that, once 
determined, topography is of no further concern as a cor- 
effect. The effective time delay measured by a 
radar system is not only influenced by the coding of the 
waveform (see, for example, Pettengill, 1971), 
but also by the radar scattering law obeyed by the surface. 
This scattering law, in general, depends on the frequency 
of the impinging radio waves, on their polarization, and on 
the angle of incidence. A generalized mapping of the 
relevant part of the surface of each target planet to the 
resolution possible with the signal-to-noise ratio available 
is an enormous, but by no means impossible, task. 

In the remainder of this subsection, we shall discuss 
the various radar methods of measurement which can yield 
information on the topography as well as the methods so far 
developed for analyzing these data to separate the (unwanted) 
planetary topographic information from the (wanted) planetary 
orbital information. 

i. Methods of Measurement 

There are a number of different methods 
that have been utilized to glean topographic information 
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from radar echoes. We will summarize very briefly three 
that relate, respectively, to topography at the subradar 
point, along the so-called Doppler equator, and over the 
visible hemisphere. The first two are being used routinely 
in planetary observations. The third has so far only been 
applied successfully to the moon since the signal-to-noise 
ratio requirements have been heretofore too severe for 
planetary application.* 

i. 1. Topography at the 

Subradar Point 

Consider a radar that 

transmits a sequence of short pulses of radio energy, or 
a phase-coded continuous-wave signal (see, for example, 

Evans and Hagfors, 1968> towards a planet. To determine 
the time delay of the reflection from the subradar point 
of these signals, a number of different techniques have been 
employed. We shall describe only one; The (decoded) echo 
power, as a two-dimensional function of time delay and 
Doppler shift, is cross-correlated with a parameterized model: 
of the expected function. As presently implemented, for 
example at the Haystack Observatory, the parameter set in- 
cludes one each for the time delay and Doppler shift corres- 
ponding to reflections from the subradar point, and others 
to describe the average, and assumed uniform, scattering law of the surface 
at the subradar point and in the surrounding regions that contribute ap- 
reciably to the echo. The values 

*A variant of this method has recently been applied with success in a radar 
experiment witli Venus the target (R.M. Goldstein, 1975) . Two antennas of 
the Goldstone Tracking Station were used to form the needed interferometer 
(see Subsection 6.i.3). ' 
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of the parameters which maximize the cross-correlation 
function are good approximations to the maximum likelihood 
estimates. The estimate for the time delay to the sub- 
radar point, of course, contains, additively, the effects 
of the topography there. Each observation, however, gives 
a delay averaged relative to that portion of the surface that 
passed through the subradar point during the observation. 

Thus, the •■footprint" on the target planet's surface for a 
delay observation depends on the duration of the observatron 
as well as on the effective pulse length and frequency 
resolution of the radar system. The minimum useful duration 
of the observation is limited primarily by signal-to-noise 
considerations, usually sufficient integration time is employed 
to insure that the interpretetion , or analysis, of the measure- 
ment will not be limited by such a consideration. For Mars, 
the integration time employed is usually far less than for 
Mercury and Venus because the rotation rate of Mars is so 
much more rapid*.. To develop this comparison quantitatively, 
consider the approximate resolution on the surface afforded 
by a pulse length (or code-element length for continuous-, 
wave signals) of and a frequency resolution of Af for 
a radar system operating at a radio frequency f . Simple 


*0f course, if sufficiently accurate knowledge of relative 
topography were available in a region, the integration could 
be extended over a longer period of time with no consequent 
sacrifice in the interpretation. 
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geometrical arguments show that 

AS = 2( 2cR (3.29) 

ir 

and 

= (3.30) 

P 

where AS represents the diameter of the spot, centered at the 
subradar point, that is "illuminated" simultaneously, or 
"encompassed" , by the pulse of length At ; and AL represents 
the length, measured along the Doppler equator*, of the 
arc whose points impart a Doppler shift to the echo that lies 
between ±Af/2 of the Doppler shift for the subradar point. 


*The Doppler equator passes through the planet's center of 
mass and is normal to the axis that results from the pro- 
jection of the "apparent" angular velocity vector onto the 
plane normal to .the earth-planet line. The "apparent" 
angular velocity vector is the angular velocity vector of 
the planet as viewed from the radar site and contains con- 
tributions from the relative orbital motions of the planet 
and site as well as from the planet's sidereal rotation. 
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The quantitites Rp and Up represent, respectively, the 
radius of the planet and the magnitude of the projection of 

X 

the apparent angular velocity vector along the Doppler 
axis. From the relation 


we can estimate the limit on integration time, t, necessary 
to prevent the surface resolution inherent in the values of 

At and Af from being "smeared" . 

Thus-, for At == 1 ysec, f - 2380 MHz, and Af - 0.1 Hz, 
which would be available initially for the upgraded Arecibo 
facility, we obtain the surface resolutions and integration- 
time limitations shown in Table 2 for each of the inner 
planets. Scaling to other values for the radar resolutions 
follows from Equations (3.29) and (3.30). Aside from the con- 
sideration of integration time, the equations show that the 
better the time resolution the smaller the area of the 
surface whose average topography affects the delay measure 
ment. The "vicious cycle" nature of this fact will be ex- 
plored in Subsection 6. ii. 

i .• 2 . Topography Along the 
Doppler Equator 

Through use of a procedure 

for data analysis different from the one described above, 
information on the relative topography along the Doppler 
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Table 2 

Surface Resolution and Integration- 
Time Limitations for Future Arecibo 
Radar Observations of the Inner Planets* 


Planet 

Surface Resolution 

Integration Time Limit 


(km) 

to Avoid "Smearing" 
(hours) 

Mercury 

77 

7.1 

Venus 

120 

16.3 

Mars 

90 

0.1 

* 

Calculation 
1 ysec and 

based on an effective 
a frequency resolution. 

pulse length. Ax, of 
Af, which matches the 


"Doppler" surface resolution, AL, to the "delay" surface 
resolution AS [see text, Equations (3.29) and (3.30) ] . 


equator can be obtained from the radar signals reflected 
from the planet. The basic idea can be outlined as fol- 
lows: The echoes, after suitable coherent integration, 

are segregated by frequency, or Doppler shift. Echoes with 
a given Doppler shift can easily be shown (see, for example, 
Evans and Hagfors, 1968) to have been reflected from a 
"strip" along the planet's surface that lies in a plane 
parallel to the one containing both the vector from the 
radar site to the planet's center and the planet's apparent 
angular velocity vector. The distance of the strip from 
the Doppler axis (which, in projection, contains the sub- 
radar point) is proportional to the relative Doppler shift 
of the strip and the subradar point. Having isolated the 
power by frequency, we may then examine it as a function of 
delay. The consequent curve will exhibit a steep rise from 
the noise level at the delay corresponding to the region near 
the Doppler equator, followed by a moderately slow decline 
in echo power for greater delays. In the idealization of 
a noise-free spherical-target situation, the first echo re- 
ceived at a given Doppler shift will be from the point on 
the Doppler equator that imparts this particular Doppler shift 
Reflections from other points on the Doppler strip will 
arrive later since they lie at greater distances from the 
radar site. If the point under discussion on the actual 
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Doppler- equator lies above or below the model spherical 
surface the echo will begin to arrive earlier or later. (A 
more accurate estimate of this time can in general be ob-. 
tained by cross correlation of the echo for the Doppler strip 
with a suitably parameterized model profile based on the 
average' scattering law for the planet.) By a comparison, 
then, of the time of arrival of the first echo from a parti- 
cular Doppler strip with the corresponding time of arrival 
of the first echo from the Doppler shift corresponding to 
the subradar point, we can deduce the relative topography of 
the two points, after correction for the spherical effect. 

In this manner, from a single observation, we can deduce the 
relative topography along an arc of the Doppler equator. 

The useful length of the arc, in general symmetrically 
placed with respect to the subradar point, will depend on 
the available signal-to-noise ratio. The scattering laws 
for the inner planets show, for example, that the echo power 
from the limb of each planet is lower than that from the sub- 
radar point by from 30 to 40 db*. Current observations have 
typically yielded useful topography for arcs of up to ±10®. 


*We exclude here observations of Venus at X-band for which 
the decline of echo power is more precipitous due to the 
added effects of the atmospheric absorption. 


-77- 


The results obtained with this approach are subject to 
distortion from several different sources: (1) An error in 

the assumed Doppler shift to the subradar point will cause 
a misidentification of the relation between Doppler shift 
and surface points along the Doppler equator with the con- 
sequence that the profile of topography will have an overall 
slope relative to the profile that should have been obtained; 

(2) The presence of anomalously high topography near, but 
not on, the Doppler equator or, conversely, the presence of 
an anomalously deep depression on the Doppler equator will 
cause the first echo for that Doppler shift to be associated 
with a point off the Doppler equator rather than on it as 
assumed in the interpretation; (3) The presence near or on 
the Doppler equator of a surface region with unusual back- 
scattering properties could lead to false 

identifications of the first echo as described above; and 
(4) Depending on the radar parameters, "aliasing", in delay 
and frequency can be present and may cause echoes from dif- 
ferent regions of the planet to be "lumped in" with the 
echoes from the strip under consideration. This fourth 
difficulty can be surmounted to a great extent by proper 
choice of radar parameters, although the size and rotation 
rate of the planet may impose certain limitations (see, 
for example, Shapiro e t al . , 1972 a, and Evans and Hagf ors , 1968) 



The other potentially distortive effects can be guarded 
against by careful comparison of the results from a number 
of observations carried out with slightly different orienta 
tions of the Doppler equator relative to the planet's 
surface and with different positions of the subradar point. 
Comparison between observations at neighboring subradar 
points can be done for Mars on the same day, because of its 
rapid rotation, as well as on different days. But for 
jijorcury and Venus whose rotation periods are comparable to 
the orbital periods, one relies on observations on adjacent 
or nearly adjacent days. Because, at a given site, obser- 
vations can extend over only a relatively small fraction of 
the day some surface regions on these two planets will pass 
through the respective subradar points when observations 
are impossible. Therefore "overlap" will be incomplete; 
nevertheless sufficient overlap exists, even in the worst 
case, to insure reliability. With the upgraded Arecibo 
facility, the useful arclengths along the Doppler equator 
will be severalfold greater and will yield more overlap . 

The comparison of the results for topography 
for points along the Doppler equator , determined from one 
observation, with the results for the same points, obtained 
as they passed through the subradar point, provides an im- 
portant check since many of the sources of systematic errors 
are different for the two methods. (Of course, these checks 
are limited to those sets of surface points that, at dif- 
ferent times,. lie both at the subradar point and at other 
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positions along a Doppler equator.) 

i. 3 . Topography Over the Visible Hemisphere 
Finally, we present a brief des- 
cription of a third technique to determine topography from 
radar observations. It is more sophisticated than either 
of the two described above, but has so far been applied use- 
fully only to observations of the moon (Shapiro et al. , 1972a; 
Zisk, 1972) . Called delay-Doppler interferometry, this method, 
at least in principle, allows the topography over a large 
area of the visible hemisphere to be determined in one 
observation, and the topography over the whole sphere 

to be determined in a series of observations. The method 
is based on the simultaneous determination of each of three 
coordinates of a reflecting region on the surface. One 
coordinate is provided by the echo delay and a second by 
the Doppler shift. These two coordinates combine to 
localize the echo to a "stick" which represents the inter- 
section of the delay contour, VA^hich is a plane normal 

to the radar site-planet vector, and the Doppler- 

shift contour,, which is a plane parallel to that formed by 
the radar site-planet vector and the planet's apparent 
angular velocity vector. The third coordinate is provided 
by interferometry. With a second antenna connected inter- 
ferometrically to the first, and receiving the same echo 
as the first antenna, one can determine the fringe phase for 



-80- 


the reflecting region isolated by the delay-Doppler co- 
ordinates. if the baseline between the two antennas, when 
projected onto the plane normal to the radar site - planet 
vector, is parallel to the delay-Doppler stick then maximum 
resolution is obtained along the stick. Quantitatively, this 
resolution is: 


AH 


RAA(J) 

2TrBp ' 


(3.32) 


where R is the earth-planet distance, X is the wavelength 

of the radio signals. Bp is the projection of the baseline 

along the delay-Doppler stick, and A<j> is the uncertainty 

in the estimate of the fringe phase. For observations of 

Venus near inferior conjunction with the upgraded Arecibo 

radar facility and its planned outrigger antenna, we would 

12 

have as typical values: R 0.3 a.u. 5 x 10 cm; 

6 

X 15 cm; B 20 km = 2 x 10 cm; and, perliaps, A(p- < 0.005 radian.* 

P 

Thus, AH - 300 m which is too crude to be very useful for 
the high accuracy required in the radar tests of relativity. 


* 

Note that the usual 2mr fringe phase ambiguity, where n 
is an unknown integer, is not a factor in these monochromatic 
observations because the a priori uncertainty in the topography 
is considerably smaller than the corresponding interval be- 
tween ambiguities. 
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If, however, the Areoibo and Haystack radars were used to 
form the interferometer and if, further, the Areoibo 
were instrumented to transmit and receive at an X-band 
frequency (X = 3 cm) , then the resolution along the Doppler 
stick afforded by the fringe-phase measurement might be 
reducible from about 300 m to the order of 10 m which would 
be useful for the relativity tests. Of course, with this 
much higher resolution the ambiguity problem* is no longer 
completely ignorable, although it is still not serious. 

The effect of fluctuations in Venus' atmosphere on the 
phase of the interferometric signals should not be serious 
(Shapiro et al. , 197^) tecause the paths from any resolution 
cell on the surface to the two radar receivers on earth 
could never be separated by more than about lo'^ radians 
and probably by one or so orders of magnitude less. Be- 
cause of the near coincidence of these paths in the atmos- 
phere of Venus, its influence on the phase delays to the 
two receivers should cancel almost completely when the signals 
are cross correlated to obtain the fringe-phase observable. 


*See footnote on p. 80. 
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The effects of the earth's atmosphere could be more 
serious; however, near cancellation can be achieved by 
use of specially chosen surface resolution cells as phase 
calibration points (Counselman et al. , 1972; Shapiro 
et al . , 1972a). 

Before concluding this discussion of the delay*Doppler- 
interferometry method of topography determination, we com- 
ment briefly on the interpretation of the three coordinates 
measurable by radar; The delay coordinate can be inter- 
preted directly in terms of a spatial coordinate; it is es- 
sentially the distance to the reflection point (for a more 
precise description, see Shapiro et al . , 1972a). The inter- 
pretation of the Doppler coordinate in terms of a spatial one 

1 

depends not only on the frequencies of the transmitted and 
received signals, but also on the (precalculated) radar site- 
planet vector and apparent angular velocity vector. The 
interpretation of the fringe-phase coordinate, as described 
above, in terms of a spatial coordinate will depend not 
only on the interpretations of the delay and Doppler coordinates 
but as well on the baseline vector that connects the inter- 
ferctneter elements and, to a lesser extent, on the (geocentric) angular 
velocity of the baseline vector. Thus only for the delay measurement is 
the interpretation in terms of a spatial coordinate reasonably direct and 
free from the need for knowledge of other .quantities. 

ii . Methods of Analysis 

We shall now discuss the 

various approaches that have been either used or proposed 
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to try to eliminate or at least to reduce the effects of topo- 
graphy on the results of the radar tests of relativity. We 
base the discussion on only the two types of delay data dis- 
cussed above: time-delay measurements to the subradar point 

and time-delay measurements corresponding to a particular 
Doppler shift different from that of the subradar point*. 

The appropriate generalization of method to include data of 


Note that this second type was considered in the earlier dis- 
cussion only insofar as it yielded topographic information 
relative to that from the subradar point. However, the second 
type of delay, like the first type, can also be considered as 
an "absolute" observable and we adopt this viewpoint here. 
There remains an important difference between the two types : 
Data of the first type can, in effect, stand alone with the 
exception noted below; data of the second type require the 
auxiliary Doppler information so that the coordinates on the 
surface to which the measurement refers can be determined. 
Aside from the consequences, discussed previously, of a mis- 
identified Doppler shift, which affects only the second type 
of delay data, both types are affected by errors in the know- 
ledge of the planet's rotation vector. At present this un- 
certainty is no problem; with improved accuracy in the surface 
resolution of the delay data, it could conceivably become a 
problem for Venus. In this eventuality the Venus radar data 
themselves could be used to improve the estimate of the 
rotation vector. 
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th© 'third "typ©, net yet available from planetary observations, 
should be obvious from context. 

The simplest approach is to ignore topography and to analyz*;; 

of the delay data based on the model of spherical planets, 
the assumption being that the topographic effects will tend 
to average out and wj 11 not seriously affect the relativity 
test results . In earlier times , when •the uncertainty in 
individual measurements of delay was comparable to the 
variation in delay imparted by topography , this approach 
was adequate. Computer experiments have shown that with the 
more accurate data already available, the correlations of 
the topography with the parameters that Characterize the 
relativity tests are high enough to have a large degrading 
influence on these tests. 

The second approach involves the use of a parameterized 
model of the topography. One can, for example, use a 
spherical harmonic expansion over the whole sphere or use 
a two-dimensional Fourier series representation confined 
to the equatorial belt on the planet that contains all pos- 
subradar points. We have tried both of these para— 
meterizations and have found the latter, augmented by 
"flattening" parameters, to be the more useful because of 
the restricted coverage of the surface afforded by the sub- 
radar points. Nonetheless, such models have a serious 
practical defect: only the relatively low frequency com— 
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ponents of the global topography can be modelled in this 
manner because of the rapid increase in the number of 
'parameters with increased resolution. Thus, to achieve 
a 10® X 10® resolution on the surface with the model we 
need, in our Fourier series approach,, about 120 parameters. 

Of course, for all three target planets combined, the 
number of topography parameters will be about 300. (This 
number can be comfortably handled in our present analysis 
program; the matrix inversion for the weighted-least- 
squares analysis for the total of over 400 parameters 
needed for the description of the orbits, etc., as well as 
of the topography, requires only a few minutes of computer 
time.) Such a low resolution, however, can not possibly ac- 
count, for example, for the myriad small craters which are 
very noticeable even with delays measured at presently achieve 
ble accuracies. To compensate for this deficiency in the 
parameterization, the errors associated with the delays in 
the analysis can be raised to a level comparable to the mag- 
nitude of the unmodelled topographic variations. Full ad- 
vantage will then not be taken of the inherent measurement 
accuracy. It appears that the radar tests of relativity 
are therefore limited by these unmodelled "high-frequency" 
terms in the topographic variations and so other techniques 
to compensate for this defect are now being applied. We 
shall describe two such techniques below. 


K 
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Cons ider first two delay measurements to a given planet 
made at widely separated times but with the same radar 
parameters and with respect to the same subradar point. From 
these two observables / two different observables can be 
formed: the difference observable and the siim observable. 

The former, of course, is simply the difference between the 
original two observables with the latter having a corresponding 
interpretation. The point of this transformation is that 
the difference observable is independent of topographical 
effects since thty enter precisely the same way in each of 
the two original observables. The difference observable could 
be used in the analysis with its full weight, considering 
only the measurement error contributions to its uncertainty. 

No allowance in the weighting need be made for the high- 
frequency components of topography being unmodelled. 

If every delay measurement had a "mate" in the sense just 
described, the topography problem would be solved in a satis- 
factory manner. The difficulty is that no delay measurement 
has a precise mate and only relatively few have near mates. 

In this situation, several types of options are open: 

(1) Ignore all observations that do not have a mate, where 
"mate" is now less precisely construed to mean an observation 
made with respect to a subradar point that falls within a 
certain specified angular separation of the subradar point 
of . the companion observation. For these "mated" observations 
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or "closure” points, analyze only the difference ob- 
servables and ignore any topographical variations that may 
exist between the two subradar points associated with the 
original observations that were used to form the difference 
observable; (2) Introduce the concept of a (possibly 
variable) correlation length for topography on each planet. 

In particular, one can continue to use a parameterized 
model for the topography, but consider all delay measure- 
ments to a given target planet to be correlated. The cor- 
relation coefficient for each such pair of delay measurements 
can be related to the angular separation of the correspond- 
ing subradar points and to either a fixed or variable cor- 
relation length. 

Both of the options described above are being inves- 
tigated. The main difficulty is the lack of a sufficiently 
large number of near mates or closure points . The basic 
geometry of radar site and planet, to say nothing of the 
difficulties in scheduling observations, severely limits 
the possibilities. For Mars the geometric cause is the 
drift in latitude of the subradar point; when observations 
are confined to short periods around oppositions, decades 
are required for all points to have near mates*. For 


With the upgraded Arecibo radar facility, observations of 
ubeful accuracy will be obtainable nearly to superior 
conjunction. But because of the limitations on the angular 
motion of the feed, observations are only possible when Mars 
is at positive declinations. 
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Mercury and Venus, slow planetary rotation and the diurnal 
motion of the earth, as well as the slow drift in latitude 
of the subradar point make it difficult to obtain closure 
points. Furthermore, the nearly eight-year periodicity in 
the relative orbital motions of the earth and Venus, and 
the nearly thirteen-year periodicity in the corresponding 
motions of the earth and Mercury, tend to introduce "aliasing" 
effects in the determination of the planetary orbits from 
the closure points. 

We must also emphasize the difficulty caused by in- 
creased measurement accuracy. Such improvements almost 
inevitably are accompcinied by increased surface resolution 
which perforce makes more difficult the obtaining of closure 
points. This vicious cycle aspect of increased measure- 
ment accuracy is partially offset by the ability to measure 
delays, in one observation, to a fairly dense set of 
points along the (instantaneous) Doppler equator. One 
might have questioned the usefulness of this technique: Why 

not simply make determinations only at the subradar point 
where the signal-to-noise ratio is highest? The answer is 

clear: Use of this "Doppler-equator" technique greatly enhances the 

number of closure points that can be obtained. 
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It is, of course, not necessary that the topography be 
determined solely from ground-based observations. Altimeters 

t 

aboard spacecraft, such as the Pioneer Venus Orbiter, can be 
used to obtain high-accuracy and high-resolution topographic 
information. These data and the radar data, where overlap 
exists, would serve to verify the accuracy of both techniques. 

Finally, to place matters in perspective, we remark that, 
at present, the radar measurement residuals attributable to 
topography, after the removal of the low-frequency terms by 
use of the two-dimensional Fourier series, range from about 
1 to 3 ysec rms for the inner planets. 

IV. Current Status and Prospects for Radar Tests of R elati\^tx 
In this final section of Part I, we give a brief sum- 
mary of the current status of the radar tests of general 
relativity and of the prospects for improvement. 
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1. Retardation of Signal Propagation 

From late in 1966 through the summer of 1971 
attempts were made using the Haystack Observatory's 
X-band (A = 3.8 cm) radar system to detect the direct ef- 
fect of solar gravity on interplanetary time delay measure- 
ments. The periods near the superior conjunctions of Mer- 
cury and Venus were especially useful for this purpose. 

Because of dwindling support for the radar system, its 
effectiveness slowly decreased with time. The combined 
analysis of these radar data yielded, through the parameter- 
ization given in Equation (2.2), a result: 

(1^) ~ 1.00 ± 0.04, (4.1) 

where the uncertainty quoted represents our best judgment > 

as to t|ie true standard error of the determination ana is 
severalfold higher than the formal standard error based on 
setting the rms of the postfit residuals to unity. The 
last published result from this series of experiments, 

[(1+Y)/2] - 1.01 ± 0.05 (Shapiro et al . , 1971b), 
was based on fewer radar data but is not substantially 
different from Equation (4.1). The main differences in the 
data sets were the additional time-delay data obtained 
near the superior conjunction of Venus in August and Septem- 
ber 1971 and the extra data obtained near the preceding 
and following inferior conjunctions which allowed improvements 
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in the determination of the topography* . 

In Figures are shown a sample of the postfit residuals for a 

typical solution involving the inner-planet radar time- 
delay data and the several hundred parameters needed to des- 
cribe the orbits, masses, and topography of these planets 
as well as various other solar“system constants. The re~ 
siduals are displayed relative to the "excess” delay given 
in Equation (2.2). Such a display is useful in that it 
shows at a' glance the relative sizes of the residuals and 
the relativistic effect, but it is also somewhat mis- 
leading in that the "masking" effects due to the correlations 
of the estimate of y with those of the other parameters is 
suppressed. One further/point must be made in connection 
with the figure and the solution on which it is based: The 

errors assumed for the measurements had a minimum value of 
3 visec, despite the fact that many of the time-delay obser- 
vations, especially near the last two inferior conjunctions 

*Because of the apparent resonance between Venus' spin and 
the relative orbits of the earth and Venus , the same longi- 
tudes , although, in general, different latitudes, are ob- 
served at inferior and superior conjunctions. 



EXCESS" DELAY {/isec) 



EARTH-VENUS TIME -DELAY 
OBSERVATIONS 

HAYSTACK | 
ARECIBO ^ 


Figure 3. Postfit residuals of Earth-Venus 
tine-delay neasurenents displaced relative 
to the excess delays predicted by general 
relativity (see text) . 
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of Venus, had measurement uncertainties as low as a few 
tenths of a microsecond. The purpose of this ad hoc 
adjustment was to compensate partially for the failure to 
model the high-frequency components of the topography 
(see Section III. 6)* complementing this radar-alone result, 
the Mariner 9 ranging data, when combined with the radar 
data yielded a preliminary value for [ (1 +t)/ 2] with only a 2% un- 
certainty. The analysis of these Mariner 9 data is being carried out cob 
laboratively between JPL and MIT; however, certain small discrepancies 
still remain to be resolved before a value can be given confidently 

and the result submitted for publication. 

In the future, at least an order of magnitude improve- 
ment in the radar-alone result for this experiment should 
be obtainable through simultaneous, or near simultaneous, 

S- and X-band observations of Mercury or Venus near 

*This particular solution did not take advantage of the 
correlations in the topography for near closure .points , 
but this neglect has no appreciable effect on either 
the estimate of y or its uncertainty: Too few closure 

pdints exist at present from the radar observations of 


Mercury and Venus, 
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superior conjunctions, using the newer, more sensitive 
radar systems. The upgraded Arecibo facility, however, will 
initially be equipped only at S-band; for this system the 
solar corona will probably limit the improvement obtain- 
able to a factor of four or so, yielding about a 2% under- 
tainty in y. If the improved antenna surface at Arecibo 
provides sufficient efficiency at X-band, then the desired 
two-freauency experiment could be conducted at Arecibo. Even 
if a high-power X-band transmitter were not available there, 
it would be possible to utilize the Haystack radar system 
for transmission of the X-band signals and Arecibo for 
i?eception, assuming, of course, that the latter can be 
equipped with a suitable receiver. Another possibility 
would be to use the 210 -foot-diameter Goldstone radar which 
is now equipped with an X-band system to complement 
its S-band one. Given the ability to carry out dual- frequency 
neasurements at Goldstone alone, or coordinated measurements at Arecibo 
(S-band) and Goldstone (X-band) , the main limitation on the accuracy 
of this experiuBnt will probably be set by the unmodelled parts of 
the topography on the target planets. Eepetition would iirprove the accuracy 

but only slowly. Nevertheless, 

there appears to us little reason to doubt that the uncer- 
tainty in y can be driven substantially below 1% through 
these radar measurements, at least down to 0.6%. 

In fact, the radar versions of this important experiment * 
■might outperform the spacecraft versions over the next 
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few years. Several circumstances could lead to this 
result: First, presently planned spacecraft involve 

ranging on two frequencies at most only on the downlink from 
the spacecraft to the ground receiver; the uplink is 
limited to S-band. Thus, depending on the spectra of the 
temporal and spatial variations of the coronal electron 
densities, it may not be possible to correct adequately for 
the coronal effects on the uplink part of the delay measure- 
ment. Second, nongravitational forces or, in the case of 
orbiters, unmodelled parts of the planet's gravity field, 
could limit the interpretation of the earth-spacecraft 
delay meaburements and, hence, the accuracy of the estimate 
of Y*- It is simply not possible on the basis of present 
knowledge to predict reliably whether the passive or active 
version of the "excess" delay experiment will yield the 
higher accuracy in the next few years. 


Landers could eliminate this problem, but the only ones 
currently planned — for Viking — have only an S-band 
capability and ranging to them will be severely limited 
by other constraints on the mission. 
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2 , Mvances of Planetary Perihelia 

The accurate determination of the advances in 
planetcxy perihelia attributable to relativistic effects 
is made difficult by the lack of an accurate, independent 
determination of the solar gravitational quadrupole moment 
as discussed in Section Ill.^.iii- 

To the extent that the parameter © that characterizes 
the solar quadrupole moment can be neglected in its in- 
fluence on planetary orbits, we can obtain a very precise 
measure of the relativistic contributions to the perihelion 
advances. As shown by Shapiro et al. (1972b), the radar 
data through 1971 implied 

= i.oOc ±0.02 (4.2) 

for the coefficient of the secular advance as given in 
Equation (2.3). The accuracy of this result relied pri- 
marily on the radar observations of Mercury. Estimates of 
a separate parameter for the secular advance of each of 
the inner planets showed that only for Mercury was the 
uncertainty under 0.1. 

With the larger set of inner-planet radar data now 
available, we have estimated the relativistic parameters 
simultaneously with J^© and the other relevant parameters 
described in the previous subsection. The preliminary 
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results indicate that 

2+2Yiij = 0.98 ±0.04 
' 3 

and 

J, = (0.5 ± 1.5) X 10"5. 

2 © 

However, because of the unusually strong correlation 
(k 0. 99 ) between these parameters, the results are very 
susceptible to significant distortion by systematic errors. 
The exposure of the effects on these results of any such 
systematic errors present in either the data or the 
theoretical model, especially of the topography, requires a 
painstaking series of sensitivity studies. In particular, 
the parameter estimates and the postfit residuals must be 
obtained for reasonably wide variations in error weightings, 
in models of the topography, etc. The results given in 
Equations (4.3) and (4.4) are based on a set of such studies 
that we regard as only partially complete; therefore no 
firm reliance can yet be placed on the numbers given; they 
are truly preliminary. 

We expect to obtain substantial improvements in these 
results not only from the continued collection of radar 
data which perforce extend the time base, but, perhaps, 
more importantly, also from the development of a much larger 
set of near- closure points. The method of determination 
of topography along an arc of the Doppler equator with each 
observation,. described in Section III. 6, has proven very 
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fruitful' due to the increase in radar system sensitivity 
recently achieved at both Arecibo and at Goldstone*. Ex- 
.amples of the topography of Venus determined by this method at 
Arecibo and at Haystack are shown in Figures 4 and 5,. 

(Campbell et al . ^ 1975) . Similar results for Mercury 
obtained at Goldstone have been published by Goldstein and 
Zohar (1974) ; however, these Goldstone delineations of the 
topography were adjusted so that the values for each 
day have a zero mean and a zero net slope. In this form they 
cannot be used to the same full advantage for the purposes 
of the radar tests of relativity as can the Arecibo and Haystack results. 

More complete coverage, with "overlap” of the arcs from dif- 
ferent days of observation, as obtained at Arecibo and Haystack, will 
enable the Goldstone data — both past and future — to be 
of full use. 

The relatively favorable 1973 opposition of Mars pro- 
vided an extended opportunity at Arecibo, Goldstone, and Haystack 

i 

to I obtain detailed and dense topographic coverage of that planet. In Figure 6 
we present a small segment of a topographic contour which 
illustrates that the terrain over certain parts of the 
surface can be extraordinarily flat. Such parts are , of 

fc 

The improvement referred to here for Arecibo relates pri- 
marily to a new feed obtained 2 years ago which yielded an 
increase in radar sensitivity of about 7 db. At Goldstone^ 
the improvement was due to the installation of a 400 kw 
transmitter on the 2 10- foot- diameter (DSS14) antenna. 
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course, highly suitable for the relativity tests since the 
correlation between near closure points can be reliably 
determined and the most effective use thereby made of the 
data for orbit analysis. Parts of the surface characterized 
by smooth, gentle undulations are almost as useful. 

We must also call attention to the order of magnitude 
difference in the uncertainties in Figure 6 as compared to 
those in Figures 4 and 5. Because of the high signal- 
to^noise ratios achievable with the Mars observations , the 
uncertainties in many of the time-delay measurements were 
at, and sometimes below, the 0.2 ysec level. With 

use of the upgraded Arecibo radar system, such accuracies 
will be achievable in time-delay measurements over major 
fractions of the geocentric orbits of Mercury, Venus, and 
Mars, and not simply in the immediate vicinities of close 
approaches. The limitation on the interpretation of the 
data for tests of relativity may then be set by the inter- 
planetary plasma unless X-band or dual -frequency observations are made 
(see Section III. 3. iii) . 

In gathering new data, attention must of course be 
directed towards securing as many "mates" as possible to form 
closure points with past data. In this manner the 

value of past data for tests of relativity can be increased 
commensurate with the intrinsic measurement accuracy 
achieved. Such attention was devoted in the 1973 Mars 
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observations and, as a resnlt, many near closure pairs 
and triplets were formed with observations from the 1971 
Mars opposition. Errors for each element of a 
pair in many cases were below 0.3 ysec. Although these 
data have not yet been fully analyzed one can safely con- 

i 

elude that the useful accuracy of these pairs for orbit 
determination will be at least at the 0.5 ysec level. This 
conclusion stands in marked contrast to the ESPo (1974) 
statement that the distance from the earth to the center 
of mass of a target planet could not be obtained by radar 
ranging to even 1 km which is equivalent to an uncertainty 
greater than 6 ysec in round-trip time delay — a limit more 
than tenfold poorer than that given above. 

To focus the above discussion on closure points and 
topography in their likely implications for improvement in 
the radar determinations of ,J and the relativistic con- 
tributions to perihelia advances, we draw the following 
conclusion: If a vigorous program for radar observations 

of the inner planets is sustained at Arecibo over the next 
lustrum, the uncertainties in the determinations of Jo and 
the combination (2+2y-3) — without the application of any a priori 

constraints based on other estimates — should be reliably bounded by 


) 




< 3 X 10 

a(2±2Yz3j , , ^ ^^-2 


(4.5) 

(4.6) 


Sufficient redundancy should be available in the coverage 
of that part of the surface spanned by the subradar points 
on each planet that topography need not be the limiting 
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factor in accuracy. The interplanetary medium would probably 
set the limit. With dual-frequency or X-band observations 
being available routinely, the limitation would be set by 
the capabilities of the timing system and of the effective pulse 
lengths planned for use initially with the upgraded Arecibo radar. Im- 
provements in these aspects, which are not yet constrained 
by the state-of-the-art, would then complete the circle by 
once again placing the accuracy limit on topography. 

We have ignored the possible contributions of Goldstone 
in the above discussion primarily because of its lack of 
availability for sustained observations, due to the demands 
for spacecraft tracking. Although it is true that the sen- 
sitivity of the Goldstone radar system will be about an 
order of magnitude or so lower than Arecibo' s, one has a 
long Way to go before system sensitivity will place the. 
primary limit on the accuracy of most of these tests of 
relativity. Therefore even limited observations at Gold- 
stone could be very useful, especially if occasionally co- 
ordinated with Arecibo 's to provide valuable independent 
checks on possible timing errors , etc.* It is rarely, 
if ever, a good policy to rely solely on one instrument for 
an important experiment. 

i Finally, we emphasize that the bounds presented in 
Equations (4.5) and (4.6) are felt to be realistic but 
not particularly conservative. Although covariance analyses 

■*The Haystack Cfcser'/atory is no longer supported for radar observations 
of the planets and so cannot contribute to such checks. 
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per formed by us indicate that substantially better ac- 
curacies can be achieved, we know from experience that un- 
accounted for systematic errors invariably degrade the 
results. 

3 . Variation of the Gravitational Constant 
The analysis primarily of five years of 
earth-Mercury radar observations yielded an upper bound 
(Shapiro et al . , 1971a) on the possible time variation of 
the gravitational constant of 
. • 

j§ < 4 X 10“^° yr”^. (4.7) 

Data accumulated over the past three years , combined with 
the earlier data, have been undergoing analyses similar to 
those described above in connection with the perihelia 
advances. In fact many of the sensitivity studies are 
identical as they serve both purposes; in others for G, 
for example, the parameters y and $ are set to unity or to 
the Dicke (1974) values 0.89 and 1.00, respectively. Pre- 
liminary results from these recent, but as yet incomplete, 
sensitivity analyses yield 

1 n“10 -1 

|g ^ ^0 yr V (4v8) 

Again we must point out that this bound is subject to the 
same caveats as were mentioned in connection with 
Equations (4.5) and (4.6). The result on the signal 
retardation, given in Equation (4 .1) , is largely exempt 
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from these difficulties for two reasons: (1) the time- 

delay effect has a unique, logarithmic signature; and 
(2) the accuracy is limited mostly by the errors in the 
measurements made near superior conjunction where the effect 


is a maximum and the signal-to-noise ratio a minimum. 

To return to the discussion of Equation (4.8), we 
note that the fourfold improvement in this bound compared 


to Equation (4.7) is due to several factors: (1) the 

increased time span of the data; (2) the improved accuracy 
of the measurements; (3) the incorporation of a model for 
the low-frequency components of the topography; and (4) the 


accurate, independent determination of the mass of Mercury. 
Until recently, the estimate of the mass of Mercury was 
dependent on the radar data -and was fairly well correlated 
(coefficient =; 0,5) with the estimate of G (Shapiro et al. , 

1971a). The Mariner 10 flyby of Mercury provided a more 
accurate, independent determination (Howard et al. , 1974): 


f 


= 6,023,600 + 600, 


(4.9) 


where M^ is the mass of Mercury in units of the sun's mass. 
This result for Mercury's mass is in embarrassingly good 
agreement with our prior publications of determinations 
based on analysis of radar data [viz. M^^ = 6,022,000 ± 53,00C 
(Ash et al. , 1967); m“^ = 6,025,000 ± 15,000 (Ash et al. , 

1971); and M^ = 6,022, 900 ± 900 (Shapiro and Reasenberg, 
1973)] and shows that our earlier result, as distinct from 
its uncertainty, was not adversely affected by the need to 
estimate simultaneously with G, 


(4.10) 


•As for the future, improvements in the estimate of 
G depend on the same factors addressed in the previous 
subsection. The same assiimptions that yielded the bounds 
given in Equations (4.5) and (4.6) also imply 

* ■ a(^) < 2 X lO"^^ yr"^ 

G 

as a reliable bound that could be achieved by the early 


198Q's, 

4. Principle of Equivalence 

There has not yet been any radar test of the Prin- 
ciple of Equivalence in regard to the relative contribution 
of the gravitational binding energy of solar-system bodies to 
their respective inertial and gravitational masses. With 
present technology, the only other possible test of this 
contribution is from laser ranging to the retrorefleotors 
•on the moon, as mentioned in' Section III.6. Radar appears 
to offer a significant test through the measurements of time 
delays of signals propagating between the earth and the other 
inner planets. In this test each of the inner planets acts 
like a small particle moving within the Jupiter-sun system, 
and, therefore, the measurements must extend in time over a 
number of orbital periods of the relevant inner planets to obtain 
an accurate result. By the early 1980-s the relative 
contribution of the gravitational binding energy should be 


determined to within a few percent if a vigorous program of 
radar measurements is pursued throughout the intervening years. 
Significantly higher accuracy, comparable to that attainable 
from lunar laser ranging, will be achievable if ranging data 
for spacecraft in planetary orbits are also obtainef. The 
best opportunities will be afforded by the Viking Mission to 
Mars and the Pioneer Venus Orbiter. Ranging data for the 
spacecraft involved in these missions could not only be of 
higher quality in their own right, but would provide unique 
opportunities to calibrate the radar measurements to these 
planets in terms of equivalent center-of-mass to center-of- 
mass delays. Independent checks of this type are very important 
■ and, in fact, enhance the value of both the spacecraft and the 

radar ranging data. 

To summarize, it is our judgment that the continuation 
of interplanetary radar measurements on a regular basis, 
especially if carried out with the equipment improvements 
suggested above for the Arecibo facility, will yield sub 
stantial improvements - from a factor of four to a factor 
of at least 10 — in each of the gravitation tests discussed in 

this section, in each case, except for the signal retar- 
dation and principle of Equivalence tests, the ground-based 
results depend vitally on observations of Mercury and there- 
fore will outperform any likely to result from spacecraft at 
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least through the end of this decade. The combination of 
the spacecraft and the radar data, of course, offers the 
highest potential . 

Further in the future one might anticipate accuracies 
sufficient to measure the next higher-order relativistic ef 
fects, but these might well require laser surveying of 
the solar system. 
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Part II 

Radio-Interferometric Measurement of the Solar Gravitational 

Deflection of Radio Waves 

I. Introduction 

. A resurgence of interest in the measurement of the 
gravitational deflection of 'light* rays by the sun has fol- 
lowed the realization (Shapiro, 1967) that radio interferonetry could 
be gainfully enployed for the purpose (Seielstad et al. , 19 70 ; 
Muhleman et al. ,1970; Sramek,1971; Hill, 1971; Sramek, 1972; 

P -.ley, 1973; Weiler et al . , 1974 ) « Here we describe in some 
detail the method of very-long-baseline interferometry (VLB!) 
that we have recently used to obtain a very accurate result 
for the deflection experiment (Counselman et al. , 1974). The 
previously published results which all involved short-baseline, 
or, more precisely, connected-element, interferometry were of 
lower accuracy;* an extended discussion of these short-baseline- 
interferometry techniques was given by Sramek (1973) . 

In our experiment, carried out in the fall of 1972 
at' a radio frequency of 8105 MHz (A 3.7 cm) , we utilized 
the 120-ft and the 60-f t-diameter antennas of the Haystack 
Observatory in Tyngsboro, Massachusetts (the "Haystack" 
and the "Westford" antennas), and two 85-f t-diameter an- 
tennas of the National Radio Astronomy Observatory (NPAO) in Green 
Bank, West Virginia, some 845 km to the Southwest. Hay- 
stack and one NRAO antenna formed a long— baseline inter- 
ferometer, and both these antennas. were directed at the 
compact extragalactic radio source 3C279. Westford and 


Recently a "long" short-baseline experiment was performed 
which yielded somewhat higher accuracy (Fomalont and 
Sramek, 1975). 
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the other NRAO antenna formed another long-baseline inter- 
ferometer and were directed at a similar source, 3C273B; 
located about 10® to the northwest of 3C279. On 8 October, 
3C279 was occulted by the sun. By observing 3C273B simul- 
taneously we were able to make relative measurements between 
the two sources , to prevent the introduction of errors by 
the separate frequency standards employed at Haystack and 
NRAO to govern the heterodyning and recording of the 
signals at the two sites. (The two antenna-receiver systems 
at a given site both utilized the same frequency standard.) 
Removal of effects of the differences between the frequency 
standards was accomplished by taking as the basic observable 
the difference between the iriterferometric fringe phases 
determined for 3C279 and 3C273B. The use of this difference 
observable also served to reduce the effects of the neutral 
atmosphere and ionosphere. 

The use of this difference fringe-phase, or DFP obser- 
vable also reduces the sensitivity to the gravitational 
deflection. The amount of the reduction is shown in Figure 7 
as a function of time. The gravitational deflection can 

I 

be detected only by determining the change in the apparent 
relative positions of the two sources during the period 



3C273 



3C279 wrt 
3C273 , 
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Figure 7 

Predicted Relative Deflections of Apparent 
Source Positions 

Numbers at various points along curves de- 
note the number of days from the occultation 
of 3C279 by the sun 


Days from 
3C279 Occultation 


A Right Ascension (milliarcseconds ) 








surrounding the occultation on 8 October. Both the right- 
ascension and the declination differences between the 
sources are changed by the gravitational deflection. 

However, because 3C273B and 3C279 both lie near the celestial 
equator r the DFP observable is sensitive mainly to their 
right- ascension difference. Iff for the purpose of explanation, 
we neglect thei solar corona, the ionosphere, the atmosphere, aber- 
ration, and numerous other small effects, we may express 
the DFP, A(}), approximately- as a function of vdie right- 
ascension difference, Aa, between the sources as 

A4) =: *Aa*cos H + A<J) (modulo 2it) (1.1) 

in which B is the length of the equatorial projection of the 
long baseline vector , A is the wavelength, and H is the hoxir 
angle at the midpoint of the baseline of a point midway between the 
sources. The quantity A4>^ is a constant which depends 
oh the baseline and the declinations of the sources but 
which, most importantly, includes the unknown constant 
difference between the phases of the independent local oscil- 
lators of the receiving stations. Because A4>^ is normally 
unknown, it is necessary to observe Acf> continuously through 

a range of hour angle H during which cos H varies signi- 

h 

ficantly, preferably near the times of rise (H a -6“) or 
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set (H » +6^), in order to determine Aa.* If the change 
of A4) between the time of rise (or set) and transit can be 
measured with an uncertainty of it follows that the un- 

certainty of the estimate of Aa is approximately 

X ^ (1.2) 

'^Aa ~ 27 tB ^(})* 

The determination of the DFP during each day of ob- 
servation without the introduction of any 2 it ambiguities, 
aside from any contained in the constant, Ac{)^ terra in 
Eq. (1.1), is a crucial requirement of the experiment. That 
is, the constant A({)^ must be the same for all measurements 
throughout a day's observations. Because the DFP is in- 
trinsically an ambiguous observable, no large gaps in its 
determination as a function of time can be tolerated. In 
particular, for a gap to be "acceptable" , .one must be 
able to connect the measurements of DFP before and after the 
gap without the Introduction of any 2ir ambiguity. During 
observations when the ray path to 3C279 passes within a 
few degrees of the sun, a gap of only a few seconds may be 


*Atmospheric phase delay increases abruptly near the times 
and another, nearly independent, determination may 

after transit on the same day. 
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unacceptable, because the solar corona introduces rapid and 
unpredictable variations in the fringe phase. On the other 
‘hand, when the separation of sun and source is 10® or 
more, even a gap of 30 min may be successfully bridged, as 

demonstrated below. 

One might wonder why the time-derivative of the dif- 
ferenced fringe phase, the so-called differenced fringe rate, 
or DFR, is not used as the basic observable since it is in- 
trinsically unambiguous. A simple calculation shows that 
the ratio of the errors in the determination of the relative 
right ascensions, Aa, for the two sources via the two methods 

is given by 

(X.3) 

Where - 27r/day is the rotational angular velocity of 
the earth and T is the integration time used to determine 
a value for the DFR. If this time is very short relative 
to a day, the accuracy achievable with the DFR observable 
may be orders of magnitude less than with the DFP, all 
other aspects of the experiment being equal. 

II. Observations 

The main characteristics of the antenna-receiver 
systems used in our experiment are described in Table 3 . 


Table 3 


Characteristics of Four-Antenna Interferometer 


Radio 

telescope 

Location 

Effective 

aperture 

(m^) 

System 

tempera- 

ture 

(K) 

Haystack* 

Tyngsboro, Mass. 

480 

80 

Westford* 

Westford, Mass. 

120 

250 

85-2t 

Green Bank, West 
Virginia 

240 

240** 

85-3t 

Green Bank, West 
Virginia 

240 

240** 

*Operated by 

the Haystack Observatory.. 


"f* 

Operated by 

the National Radio 

Astronomy 

Observatory 

single side band. 




Observed 

source 

3C279 

3C273B 

3C273B 

3C279 


(NRAO) . 
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I^cal-oscillator signals for both systems at a given end 
of the long baseline were derived from the same frequency 
standard. At the Massachusetts end, the 5 MHs signal from 
a hydrogen-maser frequency standard located at Haystack 
was multiplied to approximately 130 MHz and the resulting 
power split, one portion being sent to Westford over a 
cable whose electrical length was continuously servo- 
controlled. Although the remainder of the f requency-multr- 
plication chain at Westford was independent of that at Hay- 
stack, this additional multiplication should not have intro- 
duced a significant difference between the phases of the 
local -oscillators of the two antennas. At Green Bank, 
signal from the hydrogen-maser frequency standard was 
ried to the two antennas by buried cables. The detailed 
setup at Green Bank was essentially identical with the usual 
one employed when these two antennas form a short-baseline 
interferometer, for which the phase-stability requirements 
are far more stringent than for our experiment. The sxg- 
nals received from each source were converted from microwave 
to low frequencies, then clipped (only the sign was preserved) 
sampled, and recorded digitally at a rate of 720 kb/sec on 
magnetic tape using the HRAO Mark I system (Bare et^i^ 
i967). on each tape alternate records, of duration 0.2 sec, 
were used for a given source, so that the recorded signals 
from the two sources were interleaved. Four tape recorders 


were used at each site, in order to obtain uninterrupted 
.recording: two on line at any given time with one record- 

ing and one rewinding, and two spares. These latter were 
needed often, because the mean-time-between-f allures of 
a recorder only a few hours . Maintain- 

ing uninterrupted recordings for nine hours of observation 
each day, with each separate tape holding only three minutes 
of data, required unusual devotion of the operating person- 
nel. All told, about 5000 magnetic tapes were utilized, 
containing a total of about 3600 kilometers of tape and 

about 7 X 10^^ bits of data. 

These data were all processed at the Haystack Obser- 
• vatory with a special-purpose digital correlator connected 
to the Observatory's CDC 3300 computer. The combination 
allowed a tape pair to be processed in about four minutes, 
the output being the fringe amplitude and phase for each 
0.2-sec record; the output was then averaged coherently over 
longer intervals. Although a variety of different averaging 
intervals was used in the course of the analysis, the 
typical procedure was to form 10 -sec averages and, from 
these, the DFP observable. A simple computer program was 
used to connect the sequence of values of the DFP without the 
introduction of spurious 27T changes. However, because of 
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the subtlety of the process, we examined every single 
phase connection graphically to insure its validity; in 
every doubtful case, we re-examined the connection with 
successively smaller averaging intervals for the DFP until 
either the reliability of the connection could be assured 
or the statistical uncertainty in the determination of the 
fringe phase became too high to allow a rexiable connection. 
This latter stage was reached for an averaging interval of 
about 1 sec. In such cases we assumed the connection to be 
broken and we introduced a new unknown constant [See Eq. (1.1)1 
at the appropriate epoch into the theoretical model for the 
DFP observable. Figure 8 shows the behavior of the fringe 
phase for 3C279 for a short span of time on 4 October. 

This sample of data illustrates the obstacles to re- 
liable phase connection introduced by. data gaps and the 
solar corona. 

It is worth noting that the amount of turbulence in 
the solar corona in the relevant range of spatial scales 
(i.e., those not large compared to the Haystack-NPAO baseline) was 

exceedingly time-variable; for example, 
the DFP was very smooth for the first few hours of ob- 
servation on 3 October and then, within less than 5 minutes, 
the DFP became impossible to follow with 1-second averaging. 


^ October 1?7L 


Figure 8. Residual fringe phase as a \ 

function of time from observations of ' 

3C279 on 4 Oct 1972. These data have 
been averaged, or smoothed, over 10- 
second intervals. With this averaging 
time, the rms scatter of the phase due 
to the receiver noise and the bandwidth- 
integration-time product is approximately 20®. 
Ihe much larger phase fluctuations seen here 
are due primarily to the solar corona. 
I ! 
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Siicli ssvsire coxonsl fluctuations caused, us to eliminate 
as worthless some segments of the data. In addition, 
segments within which the phase connection was reliable, 
but which were shorter than about an hour, added no 
useful sensitivity to the determination of the relative 
position of the sources and were not utilized. 

After completion of this phase-connection and editing 
process, the DFP data were smoothed by straight-line fit- 
ting over 3-minute intervals corresponding to the duration 
of the original tape recordings to obtain one datum every 
3 minutes (except, of course, where the original data had 
contained gaps). All of the resulting data, from 23 and 
25 September, 1, 2, 3, 4, 10, 11, 13, 18, 19, and 
20 October 1972, were combined in the final analysis. 


III. Analysis 

In this analysis we assumed that the apparent solar 
gravitational deflection of the position of a source was 
given, in radians, by (Shapiro, 1967) 


GM0 

(2 + 2y)-3— / 

c P 

where Y is the unknown parameter to be estimated 
(general relativity predicts the value y = D? G is the 


(3.1) 
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gravitational constant; M is the mass of the sun; c is 

0 

the speed of light; and 

p = 2r tan (9/2) , (3.2) 

where r is the distance from the observer to the center 
of mass of the sun and 0 is the angle subtended at the ob- 
server, between the source and the center of the sun. The 
value of Y and the undeflected position of 3C279 relative 
to that of 3C273B were estimated simultaneously with a large 
set of other unknown parameters by means of iterative least- 
squares adjustment. These unknowns included the A(j>^ con- 
stants. [Eq. (1.1)] and atmospheric-model parameters, as dis- 
cussed below. 

The undeflected right ascension and declination of 
3C273B were fixed at the values given by Rogers et al. 

(1973), and the components of the Haystack-Green Bank base- 
line vector were fixed at the values determined by Hinter- 
egger et al. (1972) , with the small offsets between the 
NRAO antennas and between the Haystack and Westford an- 
tennas having been determined from a combination of conventional geodetic 
surveys and relevant short-baseline-interferometer observa- 
tions. The rotations of the baselin’e vectors with respect 
to the inertial frame of the sources were calculated using 
standard formulas for precession and nutation (Her Majesty's 


Nautical Almanac Office and the United States Naval 
Observatory/ 1961) / the instantaneous coordinates of the 

pole and UT.l interpolated between the 10-day smoothed 
values published by the Bureau International de I'Koure, 
and formulas for the displacements due to solid-earth tides 
based on Table 3a in Melchior (1966) . The 

effect of the earth's neutral atmosphere on the observable 
was calculated from the model of Chao (1968). A 
separate value of the parameter representing the zenith at- 
mospheric delay for each end of the long baseline on each 
day of observation (except, as discussed below, on 
October 4, 10, and 11) , was included as an lanknown in the 
simultaneous least-squares solution. The relatively small 
effect of the daytime ionosphere on our observable was 
calculated from a simple model in which the electron density 
varied as the cosine of the local mean solar time, and in 
which the integrated electron content along a vertical 
path at noon was taken to be 5 x 10 el m . We also model- 
led the average solar corona, assuming its density to 

vary as the inverse square of the radius from the sun, 

. -3 

With 5,000 el cm at 10 r . 

0 

The sensitivity of our gravitational-deflection result 
to the fixed value assumed for each of the parameters des- 


cribing the baseline-vector components/ the undeflected 
position of 3C273/ polar motion and the variation of UT.l, 
solid-earth tides, the ionosphere, the neutral atmosphere, 
and the solar corona was investigated by repeating the 
simultaneous solution for y, the position of 3C279, and 
all of the other unknown parameters with deliberate changes 
made to the values of these fixed parameters. We also 
examined the effects of possible 2tt phase-connection errors, 
and of omitting certain days' observations from the solution. 

IV. Result and Estimate of Uncertainty 

Our best estimate of y based on all the usable ob- 
servations, was 0.976. The- formal standard error was 
Q.009, based on the root-mean-square value of 100° of the 
post fit DFP residuals and on the assumption that the ob- 
servations, 3 minutes apart, had statistically independent 
errors. These residuals are plotted in Figure 9. It is 
clear that the errors are correlated over times significaht- 
ly longer than 3 minutes. The perhaps more accurate as- 
sumption that every 10th point, with 30-minute spacing, was 
independent would have led to a formal error greater by 
a factor of /lO, or approximately 0.03. But no estimate of 
the uncertainty which is derived entirely from the proper- 
ties of postfit residuals can be trusted: the most sig- 
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nificant, long-term trends of the measurement error as a 
function of time may have been absorbed in the model- 
fitting process in such a way that they contributed to the 
errors of the estimated parameters, but were not revealed 
in the residuals. Furthermore, errors in the assumed 
values of fixed parameters, such as the interferometer 
baseline-vector components, might seriously affect the 
solution but have no perceptible effect on the residuals. 

In order to evaluate the uncertainties due to both kinds 
of error, we performed two kinds of computer experiments 
with the DFP data. 

To estimate the uncertainty due to random, but slowly- 
varying errors, of measurement, we took the DFP data from 
each day separately, for six days (23 September, and 1, 2, 
18, 19, and 20 October) when the sun was relatively far 
from both radio sources, and we made independent solutions 
for the right-ascension difference between the sources, 
keeping both declinations fixed and setting y equal to 1. 
The rms scatter of the results of these six solutions about 
their mean value was of 00025, whereas the formal standard 
errors (assuming statistical independence for 3-minute 
spacing of the data) averaged only 0.00007. Thus, the 
results of this computer experiment confirm that the 
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originally-obtained value of the formal standard error for 
Y should be multiplied by a factor of between 3 and 4, in 
order to reflect the correlation of DFP measurement errors 
over times greater than the basic 3-minute smoothing inter- 
val. 

In order to determine the contribution to the uncer- 
tainty of our estimate of y due to errors in the values 
assumed for the parameters which were fixed in our 
solution, we repeated the original solution a number of 
times , each time changing the value of one of these para- 
meters by no less, and often by grossly more, than oiir estimate of 
its true imcertainty. The results of these sensitivity tests are 
summarized in Table 4 . A few of the entries in this table 
require explanation, as follows. 

The X and y coordinates of the pole, and the difference 
between Universal Time (UT.l) and Atomic Time (TAI , A.l, or 
UTC) are interpolated in our analysis program from values 
tabluated at 10-day intervals, with one tabular argument 
corresponding to UTC 1972 October 9.0. In the tabulation 
of each of these three quantities we made two types of 
change I (i) to the October 9.0 entry alone, and (ii) to 
the October 9.0 and all following entries equally. Thus 
we considered changes which were nearly symmetrical, and 


Table 4 . 


• Magnitudes of changes in 
by changes in the values 


the estimate of Y 
assumed for fixed 


Parameter 

Both Green Bank antennas 
geocentric radius 
longitude 

One Green Bank antenna 
radius 
longitude 
3C273B 

right ascension 
declination 

Polar motion (see text) 
variation of UT.l (see text) 

Earth tide amplitude at Haystack 
ionosphere electron density ^ 
’solar corona mean electron density 
Haystack atmosphere zenith delay 
October 4 
October 10 
October 11 


A (parameter) 

10 m 
10 m 


1 m 
1 m 

1 " 

1 " 

1 m 

2 ms 
100 % 
100 % 

100 % 

0.2 ns 
0.2 ns 
0.2 ns 


produced 

parameters 

IayJ 

0.00002 

0.00005 


0.00002 

0.00003 

0.0001 

0.00005 

0.005 

0.01 

0.007 

0.003 

0.016 

0.01 

0.006 

0.003 
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also nearly antisyinmetrical with respect to the time of the 
occultation of 3C279. In Table 4 we list the worst-case 
results, those showing the greatest effect on the estimate 
of Y. It should be noted that the estimate of y is sensi- 
tive to time-variations of the errors in the pole coordinates 
ahd UT.l, but not significantly to constant errors in these 
quantities (errors which are equivalent to constant errors 
in the baseline-vector components or the source coordinates) . 

Atmospheric-delay parameters were fixed in our solution 
only for October 4, 10, and 11 — the only days when uncer- 
tainties in the phase-connection, caused by strong solar- 
coronal fluctuations, prevented us from obtaining the un- 
broken record of the DFP over the wide range of source ele- 
vation angles which is necessary in order reliably to dis- 
tinguish atmospheric effects from source-position effects 
on the DFP. On the remaining days of our experiment, the 
phase connection was sufficiently continuous that the atmos- 
pheric parameters were able to be estimated simultaneously 
from the DFP data. The mean values of the Haystack and the 
Green Bank zenith atmospheric delay parameters obtained for 
these other days were used to fix the corresponding para- 
meter values on October 4, 10, and 11. What errors might 



-133 


have been introduced by fixing these parameters at mean 
values? our interferometric observable is sensitive prx- 
marily to the difference between the Haystack and the Green 
Bank atmospheric delay parameters. (Although both were 
normally estimated, the formal correlation between them on 
a given day was always approximately 0.99.) The rms scatter 
of the differences between the Haystack and the Green Bank 
values obtained on all the days other than October 4, 10, 
and 11 was 0.19 ns. This value provides an upper bound to 
the actual variation from day to day , because it includes 
the effects of errors in the estimates of the daily values. 
Therefore it seems conservative, in Table 4, to allow 
changes of 0.2 ns in the atmospheric-delay parameters of 

October 4, 10, and 11. 

If each of the fixed parameters listed in the first 
column of Table 4 is regarded as having a 1-a uncertainty 
equal to the corresponding entry in the second column, 
then the uncertainty in the estimate of Y due to the com- 
bination of these causes alone is given by the root sum 
qf-squares of the third-column entries, or 0.024. If, 
Xnservatively, we consider the effects of errors in these 
jiarameters to be independent of those measurement errors 
which account for the scatter of the postfit DFP residuals 


and of our six trial right^ascension solutions, we obtain 
a combined uncertainty of [(0.024)^ + (0.035) ] / = 0.042. 

Finally, we consider the possibility that, despite our 
elaborate precautions, one or more errors were committed 
in the phase-connection process, so that a spurious 2it 
change was inserted into the DFP data. For several reasons, 
the only days on which phase-connection errors could be con- 
sidered seriously are October 4 and 10 , the days 

closest to ocGultation. By deliberately in- 

serting single 2ir "errors” into the data at various times, 
we found that we could not alter the estimate of y hy as 
much as 0.01 except for October 4 and 10. We also found 
that the deletion of any single day's observations from the 
data set produced less than a 1-percent change in the es- 
timate of Y, and a negligible change in the formal standard 
error, except for October 4, 10, and (in this test) the 
llth. An accidental 2ir error on any day other than these 
three would also be rather conspicuous in the post-fit 
residuals. 

A 2ir error deliberately introduced on October 10 pro- 
duced a maximum change in y , of 0.011, if the time of the 

V. 

insertion was near 20^ UTC. On October 4 the time of 

h 

maximum sensitivity to a 2 tt error was also near 20 UTC, when 
a single 2 tt error changed the estimate of y by 0.09. For 
Such an error one hour earlier or later , the sensitivity 
was halved, and, as expected, the sensitivity approached 
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T-ero as the error- insert! on time approached the beginning or end 
of the observation span. If one assumes that exactly one 27 t 
error of either sign occurred on October 4 (on the ground that 
two or more errors of the same sign would be too conspicuous to 
go undetected), with the time of occurrence randomly distributed 
with uniform probability density between the beginning and end 
of observations, the standard deviation of the resulting error 
in the estimate of y is 0.05. 

If the observations from October 4 are deleted entirely, 
the solution for y increases by 0.069; deleting the observations 
from October 10 alone changes the solution by an approximately 
equal and opposite amount, -0.061; deletion of the October 11 
data changes the result by +0.049. The formal standard errors 
for these solutions, after multiplication by the square root 
■of ten to allow for correlations of measurement errors over 
30-minute intervals, were all about 0.05, If none of these 
days -- the only ones for which phase connection might con- 
ceivably be questionable — are included in the fit, the solu- 
tion y = 1.031 is obtained, greater by 0.055 than our all- 
inclusive result, but within one standard deviation of it. We 
conclude from these tests that coronal phase fluctuations on 

the days nearest occultation, acting both directly and in- 

1 

directly through possibly-induced phase-misconnections , 
contribute no more than 0.04 to the 1-a uncertainty of the 
combined solution for y. Combining all the sources 
of error listed in Table 4, as if independent. 
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with the coronal-fluctuation uncertainty leads to* 

Y=0.98±0.06, ^ 

. or» equivalently 

(i+1) = 0.99+0.03. 

V, Future Experiments 

Within the next few years, improved equipment and 
techniques will enable us to perform more accurate radio- 
interferometric measurements of the solar gravitational 
deflection than we were able to do in 1972. Based on the 
discussion in Section IV, we list in Table 5 all of 
the error sources which, in our 1972 experiment, we believe 
contributed more than 0.001 to the uncertainty in the 
estimate of Ih this section, we discuss, in turn, the 

most promising method of reducing the uncertainty contri- 
buted by each of these sources. 


*The results here are rounded to only two significant 
figures in view of the uncertainty which affects the 

second decimal place. 
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Table 5* Principal sources of uncertainty in the 1972 
VLB! determination of Y* 


Estimated 

source contribution to 

— Uncertainty in y 

Polar motion, UT.l, solid-earth tides 

0.013 

Atmosphere 

0.012 

Ionosphere 

0.003 

Mean solar corona 

0.016 

Inhomogeneity of solar corona, including 
possible "2 tt" phase-connection errors 

0.04 

-Scatter of postfit residuals (after al- 
lowance for correlations between errors) 

0.035 

Root-sum-of-squares of above 

0.06 
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Polar motion, variations of UT.l, and solid-earth tides 
all affect the motion of the interferometer baseline vector. 
Their effects could be determined independently, and the 


uncertainty which they introduce in the determination of 
Y reduced to negligible levels, if at least three radio 
sources were observed in addition to the one occulted by the 
sun. (Our 1972 experiment involved the observation of only 
one other source, 3C273B , by which the effects. of local- 
oscillator instabilities were eliminated.) The additional, 
baseline-calibrating, sources could be observed at night, 

before the near-sun sources had risen, or after they had 
set. 


Short-term atmospheric-delay fluctuations would not 
increase significantly if the length of the interferometer 
baseline were increased up to approximately 4000 km, for 
example by using the Haystack Observatory and the Owens 
Valley Radio Observatory (OVRO) . But the effect of these 
fluctuations in terms of equivalent source-position error 
would decrease in inverse proportion to the baseline length. 

Much beyond the 4000-km length, furbher increases would yield diminish- 
ing returns due to the decreased mutual visibility. Assuming that the at- 
mospheric fluctuations were responsible for most of the scatter of the DPP 
residuals in 1972, except on the days nearest oocultation. 
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■ we might expect a similar experiment performed with the Hay- 
stack-OVRO interferometer to yield a result for y with un- 
certainty, due to atmospheric effects alone, of about 0.007. 

Effects of atmospheric refraction could be reduced, for 
any baseline length, if additional sources which bracketed 
the position of the source occulted by the sun could also be 
observed nearly simultaneously. It is also possible that 
the atmospheric delay could be calibrated independently by 
passive radiometric measurements at infrared or microwave 
wavelengths (Schaper et aL, 1970). 

The ionosphere and the solar corona both affect our ob- 


servable by amounts which, in terms of equivalent delay, 
are inversely proportional to the square of the observing 
frequency. The most straightforward way to reduce their 
effects is to observe at a higher frequency. Several an- 
tennas are available* which are suitable for VLSI 
observations at 15 GHz, about twice the presently-used 
frequency. Connected pairs of antennas, which would enable 


the present four-antenna technique to be employed, are not 
yet equipped to operate at such high frequencies. However, 
difference observations could also be made by switching the 
pointing of a single long-baseline interferometer rapidly back 
and forth between the different sources if coronal phase fluc- 
tuations were sufficiently reduced by virtue of the higher 
operating frequency to enable phase-connection between the in- 
termittent observations of each source, or if unambiguous 


* ■ . 

These^inolude antennas at the Haystack Observatory, the NRAO 
Green Bank, the NASA Deep Space Network in Goldstone, Cali- 
fornia, and a number of others. 
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group -delay observations were made (see below) . Doubling 
the observing frequency would reduce the ionospheric con- 
tribution to the uncertainty in the estimate of y to less 
than 0.001, and the uncertainty due to the mean* solar- 
corona to about 0.004, all other aspects of the experiment 
remaining the same. 

The equivalent source-position error associated with 
a 271 phase connection error would decrease inversely if the 
baseline length were increased. Unfortunately, the proba 
bility of making 2 tt errors in the phase-connection of observa- 
tions near the sun would increase substantially , because the 
reduced correlation of coronal density fluctuations between the 
ray paths from the source to the two more widely separated 
antennas would result in greater, and more rapid, fringe 
phase fluctuations. An attempt to cope with these fluctuations 


We distinguish here the mean effect from that due to the 
inhomogeneity of the corona, which .causes short-term fluc- 
tuations in the fringe phase. 
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S imply by increasing the receiving-system bandwidth or reducing 
the system noise temperature in order to shorten the necessary 
coherent-integration time might be fruitless, for two reasons. 
First, these improvements alone would do nothing to eliminate 
interruptions of the observations due to 
equipment and operator malfunctions, occasional defects in 
magnetic tapes, etc. Rapid phase-fluctuations during such 
interruptions would cause phase-connection errors regardless 
of the bandwidth or system temperature of the interferometer. 
Second, even if all interruptions could be eliminated, and 
the phase fluctuations followed correctly, the fluctuations 
themselves might contribute approximately the same uncer- 
tainty in the determination of y a long baseline as for a 

short baseline because the size of coronal fluctuations may 
be approximately proportional to the baseline length (Knight, 1973) , all 
else remaining the same. We conclude that the use of a baseline 
much longer than that of Haystack-Green Bank might not be 
very advMitageous lonless eitl^r the observing frequency were also 



increased, to reduce the phase fluctuations due to the corona, 
or two observing frequencies were used to allow these fluc- 
tuations to be eliminated. Even with the observing 
frequency kept near the present 8-GHz value, 
however, a marginal improvement might be gained by not ob- 
serving so close to the sun. The size of the coronal fluc- 
tuations decreases approximately as the inverse square of the 
solar elongation, whereas the gravitational deflection de- 
creases only' as the inverse first power of the elongation. 

Perhaps the simplest means to surmount the problems 
of phase connection is to make use of group delay 
measurements in any future experiment. 

If the uncertainty in the determination of the group delay 
could be made less than the reciprocal of the (center) ob- 
serving frequency, then the group-delay measurement could be 
used to resolve the "2ir'' ambiguity in the associated fringe- 
phase measurement, in all observations far from the sun. 

The utilization of this "bootstrap" procedure would enable 
the full accuracy which is inherent in the use of the con- 
nected-phase observable to be achieved with only inter- 
mittent observation of a source when it- is not within a few 
degrees of the sun. When the source is close to the sun, contin- 
uous observations which enabled the fringe phase to be fol- 
lowed unambiguously by itself, and which at the same time 
included accurate group-delay measurements, would be par- 
ticularly useful because the corona introduces equal- 
magnitude but opposite-sign fluctuations into the group and 



phase delays. Thus, accurate group-delay and phase-delay data 
. can be cotibined to eliminate the effect of plasma (ionosphere as well as 
solar corona) on the observations. This iirportant possibility, if 
realized, would make it unnecessary to go to much higher 
frequencies than the presently-used 8 GHz in order to achieve 
a major improvement in the accuracy of the experiment. 

To measure the group delay at 8GHz with the requisite ac- 
curacy, error no greater than about 10 psec, will require 
substantial, but achievable, improvements in our instrumenta- 
tion. Assuming no improvement in receiving-system tem- 
peratures (~100°K) or in total spanned bandwidth of re- 
ceiver front ends (now ~ 500 MHz) , we can calculate the 
recorded-bandwidth-integration-time product necessary to 
reduce the group-delay uncertainty to 10 psec for obser- 
vations of a source with a correlated flux density of 
3 Janskys. We find, for an interferometer composed of one 
120-ft and one 85-ft -diameter antenna, that the signal- 
to-noise ratio imposes an uncertainty, a^, in the deter- 
mination of the fringe phase of approximately 

= 500 *N radians, (5.1) 

where N is the number of bits recorded at each site and 
crosscorrelated for each phase measurement. (In our 1972 
experiment, in which we recorded 7.2 x 10^ bits per second. 


this formula predicts 0.6 radians for a 1-secpnd in- 
tegration, in good agreement with our experience.) 

If we consider that the group delay will be determined from 
the fringe-phase difference between two bands 200 MHz 
apart, then the requirement of a « 10 psec leads to a 
requirement of ~ 0.009 radians for each band, or 

9 

N ~ 3 X 10 bits for each band. With our present recording 

system it would take over 2 hours to perform this 

measurement — longer than we would desire. How- 

ever, an improved recording system, called the Mark III, 
is now under development at the Haystack Observatorry and 
•promises to allow data recording rates in excess of 10 bits per 
second. Use of this system would enable the necessary 
group-delay measurement to be performed in less than 1 
minute — a very reasonable time. We note that the co- 
herent integration of the signals in each band could be carried 
out for much shorter intervals, such that solar- 
coronal fluctuations do not destroy phase coherence, and 
the between- hands phase differences from the many short in- 
tegrations could then be averaged for 1 minute or more in 
order to obtain the group delay. For example, with the 

6 

Mark III system, in 0.03 second we could record 1.5 x 10 
bits for each of two 25 MHz-wide bands, obtaining a measure- 
ment of the fringe phase in each band with « 0*4 radian, 
or of the between-bands phase difference with an uncertainty 
of 0*57 radian* The incoherent average of the 2000 
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independent measurements obtained in 1 minute would yield 
to uncertainty of 0.57/(2 000)^''^ = 0.013 radians for the 
phase difference , equivalent to 10 picoseconds in group 
delay, as required. 

we are not able to say, at present, exactly v;hat 
configuration future experiments will have, because we 
cannot predict with certainty, for example, when antenna- 
receiver systems suitable for VLBI at 15 GHz and with 
phase-stable connections between adjacent pairs of an- 
tennas, or much lower-noise or wider-bandwidth X-band front 
ends, or lO^^-bit-per- second recorders will become available 
The time scales for all of these improvements depend more on 


the availability of fmding, than 


on the electronic state 


of the art. However, it seems very reasonable to expect that, 
within 5 years, a combination of such improvements will 
enable a VLBI determination of the solar gravitational 
deflection of radio waves to be made with an uncertainty 


a (111) - 0.003 
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